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ABSTRACT 


Three  different  acoustic  phenomena  relevant  to  nuclear 
test  detection  and  diagnostics  are  discussed  in  this  report: 

(1)  The  generation  of  long  period  (3-10  minute)  acoustic-gravity 
waves  is  analyzed  in  terms  of  a  Lamb  mode  propagation  theory. 

One  result  of  the  analysis  is  that  yield-amplitude  proportionality 
is  predicted  to  break  down  for  very  large  yield  detonations. 

(2)  Tne  generation  and  propagation  of  short  period  (1-60  seconds) 
acoustic  pulses  is  treated  by  means  of  weak  shock  theory.  Yield 
and  height  of  burst  scaling  laws  are  derived  for  the  far-field 
period.  The  dependence  of  period  on  atmospheric  conditions  and 
propagation  path  is  also  discussed.  (3)  The  variation  of  long 
range  Rayleigh  wave  amplitude  with  yield  and  height  of  burst  is 
discussed  for  detonations  at  lower  altitudes  than  previously 
treated. 
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1.  INTRODUCTION 


This  report  is  concerned  with  acoustic  disturbances  produced 
by  nuclear  detonations  which  are  important  in  the  problem  of  long 
range  test  detection  and  diagnostics.  The  phenomena  we  treat,  although 
they  are  different  at  large  distances,  all  originate  in  the  detonation 
ulast  wave.  The  basic  method  of  our  treatment  has  been  to  determine 
how  the  blast  wave  characteristics,  and  hence  the  detonation  parameters 
such  as  yield  and  height  of  burst,  determine  the  far-field  acoustic 
disturbance.  We  have,  in  addition,  attempted  to  indicate  how  atmospheric 
conditions  and  propagation  path  ultimately  determine  the  far-field 
disturbance. 

Three  types  of  phenomena  are  considered: 

a.  In  Section  2  we  treat  long  period  (3-10  minute)  acoustic  distur¬ 
bances  in  terms  of  a  theory  based  on  Lamb's  atmospheric  edge  mode. 

A  detailed  consideration  of  how  the  near-field  shock  wave  evolves 
into  the  Lamb  mode  suggests  that  the  far-field  amplitude  will 
increase  somewhat  faster  than  directly  proportional  to  yield  for 
very  large  yield  explosions.  The  exact  yield  dependence  is, 
however,  determined  by  direction  of  propagation  and  atmospheric 
conditions  near  the  burst  point. 

b.  In  Section  3  we  discuss  the  dependence  of  short  period  (1-60  second) 
acoustic  signals  on  source  and  propagation  parameters.  We  find 
that  the  signal  period  varies  approximately  as  the  cube  root  of 
yield,  is  relatively  independent  of  burst  height,  and  is  very  much 
affected  by  the  direction  of  propagation  relative  to  the  high 
altitude  winds. 

c.  In  Section  4  we  consider  the  variation  of  far-field  Rayleigh  wave 
amplitude  with  yield  and  height  of  burst  for  atmospheric  explosions. 
This  section  is  primarily  an  extension  of  previous  work  to  lower 


heights  of  burst  than  were  formerly  considered.  Our  principal  new 
finding  is  that  the  amplitude  for  low  altitude  detonations  is  a  much 
more  rapidly  increasing  function  of  yield  than  for  higher  altitude 
detonations. 

Overall  conclusions  reached  are  presented  in  Section  5. 
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2.  SOURCE  MODELS  FOR  THE  EXCITATION  OF  LAMB'S 
ATMOSPHERIC  EDGE  MODE  BY  NUCLEAR  EXPLOSIONS 


2.1  Introduction 

The  acoustic  waveform  observed  several  thousand  kilometers 

away  from  a  megaton  range  nuclear  detonation  begins  with  several 

cycles  having  a  dominant  period  of  from  about  3  to  10  minutes.  The 

amplitude  of  this  portion  of  the  acoustic  signal  has  previously  been 

predicted  to  be  proportional  to  explosion  yield^  Recently  these  first 

(2) 

few  cycles  of  the  waveform  have  been  analyzed  by  Pierce  and  Posey 

under  the  hypothesis  that  the  acoustic  pulse  propagates  as  the  real 

atmosphere's  counterpart  of  Lamb's  edge  mode.  In  the  material  which 

follows  we  rely  heavily  on  the  Pierce-PoSey  theory.  This  theory 

(3) 

appears  to  be  in  at  least  qualitative  agreement  with  test  data. 

When  it  has  been  tested  against  multimode  numerical  calculations  the 

theory  has  also  produced  good  agreement  for  the  first  few  cycles  of 

(2) 

the  waveform.  For  our  purposes  the  significant  feature  of  this 

theory  is  that  it  is  analytic.  This  allows  us,  for  the  first  time, 

to  analyze  in  some  detail  the  relationship  between  the  near-field 

shock  wave  and  the  resulting  far-field  waveform. 

Specifically  we  ask  the  following  question:  At  what  distance 

from  the  detonation  and  in  what  direction  should  one  utilize  the  shock 

wave  parameters  to  begin  a  linear  propagation  theory?  Obviously  a 

linear  theory  is  inadequate  viien  the  shock  is  strong.  For  a  1  MT 

detonation  at  sea  level,  for  example,  the  shock  becomes  weak  (the 

(4) 

relative  overpressure  is  .1)  at  a  distance  of  about  9  km.  Since 
this  is  larger  than  the  sea-level  atmospheric  scale  height,  the  effects 
of  atmospheric  inhomogeneity  might  reasonably  be  expected  to  play  a 
role.  In  other  words,  the  weak  shock  parameters  such  as  relative 
overpressure  and  positive  phase  duration  will  be  different  in  different 
directions.  Therefore  we  want  to  know  whether  a  particular  portion  of 
the  shock  front  dominates  excitation  of  the  Lamb  mode  and,  if  so,  we 
want  to  know  what  the  shock  parameters  are  for  this  portion  of  the  front. 
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The  question  of  what  parts  of  the  shock  front  are  most  impor¬ 
tant  in  determining  the  far  field  waveform  is  rot  the  only  question  we 
have  to  answer.  The  distance  at  which  we  match  the  weak  shock  to  a 
linear  propagation  theory  will  also  make  a  difference  in  the  far-field 
result  which  the  theory  predicts.  For  example,  consider  the  differ¬ 
ence  that  occurs  if  we  take  the  horizontally  propagating  shock  wave 

from  a  1  MT  sea-level  explosion  at  two  places:  (1)  9  km,  where  the 

(A) 

relative  overpressure  is  .1,  and  (2)  58  km  where  it  is  .01.  For  a 
weak  shock  in  a  homogeneous  atmosphere  irRt+  is  a  constant. n  Is 
the  peak  relative  overpressure,  R  is  the  distance,  and  t+  is  the  posi¬ 
tive  phase  duration.  In  the  present  example  t+  must  increase  by  a 
factor  of  about  1.55  between  9  km  and  58  km.  The  energy  in  the  low 
frequency  portion  of  the  pulse  >>  t+j ,  which  is  responsible  for 
the  long  period  portion  of  the  far-field  waveform,  is  proportional  to 
ifRt  2.  This  result  follows  from  Fourier  analysis  of  a  weak  shock 
profile  such  as  the  so-called  Glasstone  pulse.  Since  nRt+  is  a 
constant,  the  low  frequency  energy  is  proportional  to  t+  and  must  also 
increase  by  a  factor  of  1.55  between  9  km  and  58  km.  We  conclude  that 
if  we  chose  to  match  a  linear  propagation  theory  to  the  shock  wave  at 
n  *  .01  rather  than  at  .1,  we  would  predict  far  field  amplitudes  which 
are  greater  by  a  factor  of  1.55  The  difficulty  is  not  that  wc  have  not 
gone  to  small  enough  tt  .  For  a  weak  shock  n  never  becomes  so  small 
that  the  propagation  is  linear;  that  is,  the  Fourier  amplitudes  at 
different  frequencies  are  never  frozen  relative  to  each  other.  The 
energy  in  the  low  frequency  portion  of  the  pulse  continues  to  increase 

proportional  to  t  which  quantity  increases  asymptotically  as  *4nK 

+  (5V 

for  a  horizontally  propagating  shock.  ' 

In  the  sequel  we  use  the  Pierce-Posey  theory  to  investigate 
how  the  far-field  waveform  evolves  from  the  near-field  shock.  Specific 
conclusions  reached  a.. a:  (a)  The  shock  front  extending  over  a  number 
of  scale  heights  in  altitude  is  of  almost  equal  importance  in  exciting 
Lamb's  atmospheric  edge  mode.  A  consequence  of  this  fact  is  that 


features  peculiar  to  a  limited  portion  of  the  shock  front  cannot  domin¬ 
ate  the  excitation.  An  example  of  what  we  mean  by  a  feature  peculiar 
to  a  limited  portion  of  the  shock  front  is  the  existence  of  a  shock 
precursor  region  near  the  ground  within  which  shock  parameters  are 
drastically  altered.  This  occurs  when  a  low  altitude  detonation  occurs 

over  a  heat  absorbing  surfaced  (b)  For  small  yields,  the  appropriate 

1/3 

distance  for  matching  increases  as  Y  ,  where  Y  is  the  explosion 
yield.  For  large  yields  the  appropriate  distance  for  matching 
increases  more  rapidly  than  Y^^.  This  produces  a  yield  dependence 
in  the  far-fleld  amplitudes  which  varies  approximately  as 

f  _  - - 

V  |  1  +  ~  "V  n  7*  j  for  Y  >  Y*,  and  as  Y  for  Y  <  Y*.  Typically,  Y*. 

which  depends  on  atmospheric  conditions  near  the  burst  nnd  on  the 
direction  of  propagation,  is  the  order  of  10  Ht. 

In  the  next  section  we  outline  the  method  used  and  the  maj  or 
assumptions  made  in  calculating  the  far-field  Lamb  mode  amplitude  from 
the  near  field  stuck  parameters.  The  calculation  itself  is  contained 
in  an  appendix. 

Then  in  Section  2,3  we  compare  the  above  theoretical  prediction  (b) 
with  the  test  data.  While  the  results  of  this  comparison  are  not  con¬ 
clusive  we  find  that  there  is  indeed  a  strong  indication  that  on 
occasion  very  large  yields  do  produce  anomalously  large  far-field 
amplitudes. 

Finally,  we  summarize  the  results  of  this  study,  and  indicate 
where  we  believe  extensions  of  the  analysis  would  be  most  fruitful,  in 
the  conclusions  presented  in  Section  2,4. 

2 . 2  Method  of  Calculation  of  the  Far-Field  Lamb  Mode  Amplitude 

In  this  section  we  outline  the  method  used  and  the  major 
assumptions  made  in  calculating  the  fat  field  Lamb  mode  •mplitude.  The 
actual  calculation  is  done  in  Appendix  A. 

The  fundamental  assumption  of  the  theory,  the  justification 
for  which  has  been  discussed  by  Pierce  and  Posey,  is  that  the  long 
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period  acoustic  pulse  observed  at  great  distances  propagates  in  the 
real  atmosphere's  counterpart  of  Lamb's  edge  mode. 

To  extract  this  mode  from  the  near  field  shock  wave  we  pick 
a  cylindrical  matching  radius  rQ,  in  a  manner  described  below,  and 
decompose  the  shock  pressure  pulse  into  a  modal  sum,  one  of  whose 
terms  is  the  desired  Lamb  mode.  This  Lamb  mode  term  is  then  taken  to 
propagate  linearly  to  the  observation  point,  a  great  circle  distance  r 
from  the  burst,  in  accordance  with  the  Lamb  mode  dispersion  equation. 
Because  we  arc  primarily  interested  in  modeling  the  source,  we  do  not 
attempt  to  incorporate  a  realistic  atmospheric  model  in  the  theory, 
i.e.,  one  which  wmld  include  effects  such  as  acoustic  ray  refraction. 
Horizontal  refraction  of  the  Lamb  mode  may  be  an  important  effect,  but 
this  can  be  treated  independently  from  the  source  modeling. 

To  simplify  the  analysis  several  assumptions  arc  made  con¬ 
cerning  the  properties  of  the  near  field  shock  wave: 

a.  The  weak  shock  is  characterized  as  a  Glasstone  pulse  wliosc  peak 
overpressure  s  and  positive  phase  duration  t+  arc  taken  to  be 
functions  of  altitude  and  range  as  described  below.  It  must  be 
admitted  that  the  choice  of  a  Glasstone  pulse  is  made  for  lack  of 
any  better  analytic  representation  of  the  detonation  shock  wave. 

To  our  knowledge  no  detailed  analysis  lias  ever  been  made  of  how 
adequately  the  Glasstone  pulse  represents  the  very  low  frequency 
components  of  the  detonation  shock  wave. 

b.  Based  on  test  observations  it  is  assined  that  t+  is  much  smaller 
than  the  far-flold  period  T.  Use  is  made  of  thin  assumption  in 

expanding  the  far-field  amplitude  as  a  scries  in  ~  and  keeping 
only  the  first  non-zero  term. 

c.  The  shock  front,  except  where  it  intersects  the  ground,  is  assined 
to  be  spherical.  This  assumption  is  only  intended  to  apply  to  the 
portion  of  the  shock  which  Is  effective  in  exciting  the  Lamb  mode, 
it  need  not  apply,  for  example,  to  the  upward  going  shock  front. 
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This  assuaption  will  be  correct  if  the  shock  becomes  essentially 
sonic  within  an  effective  scale  height  in  the  direction  in  question. 

d.  For  low  altitude  detonations  the  secondary  shock  which  results  when 
the  direct  shock  is  reflected  from  the  ground  is  assumed  to  make 
a  contribution  to  the  far-field  amplitude  equal  to  that  of  the 
direct  shock.  The  justification  for  this  assiraption  is  that  for  a 
low  altitude  detonation  the  secondary  shock  at  sufficient!)  large 
distances  from  the  burst  point  has  essentially  the  same  amplitude 
and  radius  of  curvature  as  the  direct  r.hock.  Furthermore,  the 
secondary  sliock  follows  the  direct  shock  within  a  time  which  is 
small  compared  to  the  periods  observed  in  the  far  field.  For  a 
contact  or  surface  burst  the  secondary  and  direct  shock  arc  of 
course  the  same,  although  the  amplitude  may  be  reduced  due  to 
energy  expended  in  digging  a  crater,  etc.  This  assunption  restricts 
the  present  theory  to  low  altitude  detonations  (less  than  say  IS  mt). 
For  sufficiently  high  altitudes  the  radius  of  curvature  of  the  two 
shocks  may  be  very  different  and  the  time  lag  between  the  shocks 
may  be  comparable  to  the  far  field  period.  In  this  case  we  would 
expect  interference  effects  to  occur.  The  high  altitude  case  has 
not  teen  investigated. 

The  result  of  the  preceding  steps  and  assumptions  is  Eq.  (A-30) 

of  the  Appendix: 


r 


L 


Y 


«*> 


/ 


o 


r 

V>*2  ♦  r2  -  r  1 

V  o  o 

Z 

• 

-t  + 

e  rH  *(r  ,z)t2  (r  ,z)  dz 

0  ▼  o 

(A-30) 

V* 

c  'a  J 

• 

-  7  - 


where: 


■  pressure  amplitude  of  the  far  field  Lamb  mode 

P  •  sea  level  pressure 

s 

tq  ■  cylindrical  matching  radius  (to  be  discussed  below) 

r  ■  earth  radius 
c 

r  ■  great  circle  distance 

Y  ■  ratio  of  specific  heats 

H  ■  scale  height 

t  .  -  characteristic  dispersion  time  which  is  approximately  a  quarter 

of  the  far  field  period  T. 

A ^  •  derivative  of  Airy  function  with  respect  to  its  argument, 
t  -  time  after  arrival  of  the  far-field  acoustic  signal 
c  ■  sound  speed 

Az  ■  vertical  distance  relative  to  height  of  burst 
z  ■  altitude  relative  to  sea  level 
H  -  peak  overpressure  of  weak  shock 

t+  ■  positive  phase  duration  of  weak  shod. 

To  proceed  from  Lq.  (A-30)  it  is  necessary  to  assume  a  specific 
spatial  dependence  for  “  and  t+.  We  take  t+  to  be  given  at  a  small 
radius  Rq,  where  n (Rq)  ■  .1,  by  the  problem  M  results  for  a  homogeneous 
atmosphere  with  ambient  pressure  equal  to  the  ambient  pressure  at  the 
burst  point.  Rq  is  also  taken  from  the  problem  M  results  and  both 
t+(RQ)  and  Rq  are  scaled  according  to  Sachs  scaling.  At  the  spherical 

I  2  2 

radii  /Az  +  ro  >  Rq  where  t  and  t+  arc  required  in  Eq.  (A-30)  these 

quantities  arc  obtained  from  the  values  at  K  by  using  Reed's  theory 

°  (5) 

for  weak  shock  propagation  in  an  exponential  atmosphere.  After 
some  labor  it  is  found  that  Eq.  (A-30)  can  be  reduced  to  be  written 
solely  in  tenus  of  the  horizontal  values  of  the  shock  parameters  *  aid 
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(A-44) 


(r  ,z)  r 
o  o 


AJ 


,1 


wlierc  /  _Is  the  height  of  burst  and  the  new  integration  variable 

-v  /.  2  ^  2 

\  Az  +  r  -  r 
V  o  o 

y  -  - 


CT 


i quat ion  (A-44)  is  essentially  identical  to  the  result  of 
Pierce  and  Posey.  This  is  in  spite  of  the  fact  that  we  have  used  a 
distinctly  different  shock  nodel,  one  which  includes  the  effects  of  an 
exponential  atmosphere.  This  result  occurs  because  in  the  Integral  of 
Lq.  (A-30)  the  effect  of  the  exponential  atmosphere  on  the  shock 
parameters  tends  to  be  cancelled  by  tho  Lamb  mode  weighting  factor 
e  Z^Y**  as  well  as  other  purely  geometrical  effects. 

Now,  the  first  maximum  of  the  far-fleld  amplitude  | —  j 

a;  C.-iarJ 

occurs  at  about  t  ■  t^.  For  t  •  the  quantity  occurring  in 

Eq.  (A- 39)  is  shown  to  be  significant  only  between  about  0  <  y  <  2. 
Truncating  the  y  integration  at  y»2  corresponds  to  an  effective  upper 
limit  on  the  Az  values  of  interest  of  Az  £  20^,  that  is  about  a  half 
wavelength  of  the  far-field  disturbance.  Typically  this  means  that 
the  shock  at  Az  values  of  at  least  up  to  60  km  contribute  almost  equally 
to  the  Lamb  mode  excitation.  This  is  the  basis  for  our  statement  that 
features  peculiar  to  a  limited  portion  of  the  shock  front,  such  as  the 
existence  of  a  precursor  region  near  the  ground,  cannot  dominate  the 
Lamb  mode  excitation. 

It  now  remains  to  choose  the  matching  radius  tq  which  occurs 
in  Lq.  (A-44) .  Physical  argwents  are  given  in  the  Appendix  that  this 
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should  bo  done  in  the  following  way:  Ue  make  a  Fourier  analysis  of  the 
Glasstonc  pulse  and  assume  that  each  frequency  component  within  tin 
pulse  propagates  at  »  group  velocity  given  by  the  Lnmb  mode  dispersion 
equation.  Tltc  matching  radius  rQ  is  determined  by  requiring  that  at 


the  low  frequencies 


u  a*  which  will  be  important  in  the  far 

d 


field  must  Itave  propagated  ahead  of  the  largest  amplitude  frequency 
components  of  the  weak  shock  (u  ^  * 


w 


The  basis  for  this 


criterion  is  that  while  the  low  frequency  components  are  located 
behind  the  shock  front  they  will  gain  energy  from  the  higher  fre^  icy 
components,  but  once  they  move  ahead  of  the  higher  frequencies  their 
propagation  should  be  amenable  to  a  linear  treatment.  A  more  detailed 
statement  of  this  argument  is  given  in  the  Appendix. 

This  criterion  assises  the  form: 


t+(ro)  -  td(ro) 


(A-54) 


which  serves  to  determine  rQ.  liquation  (A-54)  has  an  additional  inter¬ 
pretation;  Suppose  we  start  with  the  near  field  disturbance  and  follow 
its  development  to  increasing  radii  by  means  of  weak  shock  theory.  In 
addition  suppose  we  take  the  far-fleld  disturbance  and  carry  it  back 
toward  the  detonation  by  means  of  the  Lamb  mode  theory.  According  to 
Eq.  (A-54)  the  matching  radius  is  the  distance  at  which  the  dominant 
frequency  content  of  the  two  disturbances  becomes  equal. 

It  turns  out  that  for  small  yields  t^fr)  >  t+(r)  for  nil  r 
at  which  the  shock  is  weak.  In  this  case  we  use  conventional  Sachs 


scaling  r  a 
o 


Y_l 

Fb! 


1/3 


,  where  Y  Is  the  detonation  yield.  The  constant 


of  proportionality  is  chosen  to  make  r 


Y 

T 


a  continuous  function  in 


making  the  transition  to  the  large  yield  case  where  hq.  (A-54),  rather 

than  Sachs  scaling,  is  used  to  determine  r  . 

o 

The  detonalioa  shock  parameters  enter  into  Eq.  (A-44)  in  the 
-  -  2 

form  n(r^,z)t+(r  ,z)r  .  Since  for  a  weak  shock  the  quantity 
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.  _  IV  \2/3 

a(r,z)t+(r,z)r  is  proportional  to  J~*|  and  independent  of  r  this 

means  that  the  far-field  amplitude  is  just  proportional  to 

•  y  2/3 

—  t.(r  )  which  increases  «ore  rapidly  than  when  r  is  determined 
1 b  °  Pb  ° 

by  tq.  (A>*54).  The  result  of  this  analysis  is  the  following*: 
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(A-69) 


where  P.  is  in  atmospheres,  Y  in  kT,  and  cr  .  and  r  in  km.  The 
l  a  c 

"correction  factor"  Q  is: 


qJL-U  ,  H 

Q;pb  3  Y*  pb 


A 


Y__  Y* 

*  pb  -  pe 


1  *_  <  X!!. 

pb-pb 


(A-75) 


The  critical  yielJ  to  pressure  ratio  is  given  in  terms  of  Pierce  and 
Posey's  dispersion  parameter  h^  by: 


sec 


3/2 


Xi  Aorwh  Mt _ 

P*  hkk^  atmosphere  ^9/2 


(A-71) 


*  The  corresponding  constant  in  th»  work  of  Pierce  and  Posey^2)  would 
appear  to  be  3.2  x  10“3  rather  than  our  value  of  4.7  x  10“3.  This 
difference  probably  arises  primarily  from  their  choice  of  Wt?  r0  as 
(3.4  x  10*2) (.33)2(1.61)  •  6.1  x  10“3  km/sec2  (or  a  1  kT  sea  level 
detonation.  Based  on  problem  M^)  we  lave  taken  nt2ro  ■  ( .l)(l/3)2(  ,9)*10 
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The  qu  int  il*  li  is  i  dependent  of  fr«*';uenc\  hut  Is  a  function 

of  atmospheric  conditions  nyr  the  Ixirst  (we  only  ropiirc  li  in 

l.q.  (A-hh)  out  to  tin*  distance  r  )  as  well  as  tin*  direction  of  |  ropa- 

°  v* 

ration.  For  this  reason  it  appears  that  an  only  he  calculated 

1  I 

numerically  in  rnnjun<  lion  with  a  detailed  atmospheric  model  for  enndi- 
t  ions  near  the  hurst  point  . 

Nevertheless,  as  discussed  in  the  next  section,  it  is 
possible  to  set.  some  evidence  in  the  low  alliluci  lest  data  of  the 
breakdown  of  yicl  d- impl  itude  proportionality  for  very  large  yield 
del  tinat  ions. 


2.  i  Com  par  Ison  of  Tuci  r  with  Test  data  i  ir  Cargo  Yields 


There  Is  at  least  a  suggestion  in  the  lost  data  that  very 
large  yields  do  on  occasion  produce  anomalously  large  far-l icld  ampli- 

I’j  (t“  ([)  -\J  sin  r/r^  ^ 

tudes.  According  to  la|.  (A-64)  the  quantity  — 1 - *  y""o(Y) - - 

is  sup|xi.scd  to  lie  i  constant  for  near  surface  delon.it  ions.  l  aking 
l'j(t»ij)  to  In  ahuit  Ikilt  i  he  pressure  amplitude  lirst  peak  to  through 
Pf  ,  flnd  Tj  to  ,m;  about  a  quarter  of  the  far-field  period  T,  we  have 

I’,  t  'J  sin  r/r  j/ » 

used  the  data  presented  in  Ref.  4  to  plot  Y  vs.  — —  1  t 


in  Fig.  2.1,  Note  lint  the  plot  is  somi  logarithmic,  ilie  detonations 
used  are  the  Soviet  explosions  of  (a)  10  September,  (h)  11  September, 
(c)  14  September,  (d)  «  October,  (e)  6  October,  (f)  20  October, 

(g)  23  Octobir,  (h)  10  October,  (i)  31  October,  I9ftl,  and  the  U.  S. 
explosions  of  (j)  4  May,  (k)  10  June,  (1)  12  June,  (m)  27  June,  and 
(n)  II  July,  1901.  Some  detonations  occur  In  Fig.  2.1  more  than  once 
because  data  fron  several  stations  are  used. 

The  vertical  lines  in  Fig.  2.1  have  been  arbitrarily  drawn  at 
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FIGURE  2-1  s  OBSFRVEn  DEVIATION  FROM 

YIELD-AMPLITUDE  PROPORTIONALITY 
FOR  LARGE  YIELDS 


I 


y  r  t3/2  |h^t  .  (100  sec)3/2 

_ I _ J 

AO  (To 


l3 


bdi*  3/2 

of  25  M^-ar  (100  sec)  corresponds  to  a  choice  of  sound  speed 
c  =  .  31  km/sec  in  Lq.  (A-64) .  The  detonations  which  lie  outside  the 
large  amplitude  boundary  do  tend  to  be  large  yield  detonations.  For 
example,  58  and  24  Mt  at  two  stations  lie  outside,  25  Mt  and  two  9Mt 
detonations  also  lie  outside. 

Figure  2.1  provides  some  indication  of  the  possible  breakdown 
of  yield -amplitude  proportionality.  However,  we  should  caution  that 
some  of  the  deviations  may  be  caused  by  propagation  rather  than  source 
effects  and  furthermore  that  some  of  the  yield  estimates,  which  are 
based  on  seismic  data,  may  be  in  error. 


2.4  Conclusions 


We  have  found  that  the  shock  front  extending  from  the  ground 
to  many  scale  heights  above  the  ground  is  of  almost  equal  importance 
in  determining  the  far-field  Lamb  mode  amplitude.  Because  of  this, 
the  far-field  amplitude  should  be  relatively  independent  of  phenomena 
peculiar  to  any  small  portion  of  the  shock  front.  In  particular, 
ground  effects  such  as  the  presence  of  a  precursor  should  not  play 
a  very  important  role. 

For  large  yield,  low  altitude  detonations  we  have  found  a 


far-field  amplitude  dependence  which  varies  as  Y 


where  typically  Y*  is  the  order  of  10  Mt  for  a  sea  level  detonation. 
This  dependence  arises  from  a  detailed  consideration  of  the  mechanism 
by  which  the  low  frequency  components  cf  the  near-field  (shock) 
disturbance  become  spatially  separated  from  the  high  frequency  com¬ 
ponents.  The  precise  value  of  the  critical  yield  Y*  where  yield- 
amplitude  proportionality  breaks  down  is  determined  by  atmospheric 
conditions  near  the  burst  as  well  as  the  direction  of  propagation  and 
must  be  determined  by  detailed  numerical  calculations.  Nevertheless, 
in  analyzing  the  test  data  we  have  seen  some  evidence  that  very  large 
yields  do  on  occasion  produce  anomalously  large  far-field  amplitudes. 
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Finally,  we  note  that  a  number  of  potentially  important 
source  effects  can  be  explored  within  the  framework  of  the  present 
theory.  The  most  important  of  these  effects,  which  we  have  not  investi¬ 
gated,  are  probably  related  to  the  presence  of  a  secondary  shock  front 
due  to  ground  reflection.  We  expect  that  at  low  altitudes  the  effect 
of  the  secondary  shock  in  determining  the  far-field  amplitudes  is 
approximately  equal  to  that  of  the  direct  shock,  but  for  higher  alti¬ 
tude  detonations  interference  or  destructive  effects  between  the  two 
shocks  may  occur.  It  would  be  desirable  to  subject  these  ideas  to  a 
quantitative  analysis. 
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3.  VARIATION  OF  FAR-FIELD  HIGH  FREQUENCY  ACOUSTIC 
PERIODS  WIT!!  YIELD  AND  HEIGHT  OF  BURST 


3.1  Introduction 

This  portion  of  the  report  is  concerned  with  detonation- 
produced  infrasonic  signals  which  have  periods  between  about  1  second 
and  1  minute.  These  signals  have  been  far  less  extensively  studied 
than  the  long  period  (3-10  minute)  acoustic-gravity  waves  discussed 
in  Section  2. 

The  short  period  disturbances  are  difficult  to  treat 
theoretically  because  their  propagation  is  sensitive  to  small-scale 
meteorological  phenomena.  For  example,  Meecham  has  calculated 

that  signals  in  the  second  to  minute  period  range  have  a  substantial 
probability  of  encountering  a  diffracting  wind  duct  in  propagating 
between  the  ground  and  an  upper  boundary  at  about  50  km  altitude.  On 
each  encounter  the  acoustic  signal  is  split  and  a  small  fraction  of 
the  energy  is  lost  from  the  main  pulse.  Accord. nc  to  Meecham,  since 
the  propagation  path  to  large  distances  involves  a  large  number  of 
transits  between  the  ground  and  the  upper  boundary,  the  main  pulse 
eventually  becomes  so  degraded  at  large  ranges  that  it  becomes  lost 
in  the  multiplicity  of  pulses.  In  fact,  as  Meecham  notes,  the  far- 
field  signal  is  observed  to  consist  of  a  large  number  of  apparent 
pulses  extending  over  a  time  interval  of  an  hour  o’-  more.  He  attri¬ 
butes  the  extended  signal  duration  to  horizontal  refraction  by  large- 
scale  weather  fronts. 

Meecham' s  findings  are  related  to  the  degradation  of  pulse 
amplitude  with  range.  We  are  concerned  with  the  far-field  pulse  period 
and  here  the  situation  is  quite  different.  We  believe  that  the  period 
of  the  main  pulse  is  relatively  unaffected  by  pulse  splitting  phenomena. 
In  particular  we  find,  as  outlined  in  the  subsequent  sections,  that 
the  disturbance  period  is  essentially  determined  after  only  a  few 
transits  between  the  ground  and  the  upper  boundary;  that  is,  before 
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pulse  splitting  has  had  an  appreciable  effect.  A  corollary  to  this 
finding  is  that  the  far-fleld  periods  have  a  very  weak  dependence  on 
range. 

The  model  on  which  the  analysis  is  based  follows  Meechaai  in 

describing  the  upper  boundary  as  a  reflecting  plane  at  about  50  km 

altitude.  This  is  the  mode  of  long-range  propagation  observed  for 

(9) 

infrasound  from  rockets.  The  precise  altitude  of  reflection,  which 
is  an  important  parameter  in  the  theory,  is  determined  prisarily  by 
the  direction  of  signal  propagation  relative  to  the  winds  at  about 
50  km  altitude. 

in  our  model  the  period  observed  in  the  far-fleld  is  pro¬ 
portional  to  the  positive  phase  duration  of  the  near-field  shock  wave. 
The  positive  phase  duration  increases  as  the  shock  propagates  and  this 
effect  is  calculated  by  means  of  weak  shock  theory.  This  is  the  novel 
feature  of  our  analysis,  namely  that  it  includes  nonlinear  weak  shock 
effects  over  large  distances. 

To  obtain  starting  values  for  the  application  of  weak  shock 
theory  we  use  modified  Sachs  scaling.  This  procedure  has  been  demon¬ 
strated  to  be  superior  to  ordinary  Sachs  scaling  for  prediction  of  peak 
overpressures  in  an  inhomogeneous  atmosphere.  We  extend  the  concept 
of  modified  Sachs  scaling  to  apply  to  prediction  of  positive  phase 
duration  as  well  as  peak  overpressure. 

In  Section  3.2  we  develop  the  weak  shock  propagation  model 
and  in  Section  3.3  we  obtain  starting  values  for  this  model  using 
modified  Sachs  scaling.  The  theory  is  then  used  in  Section  3.4  to 
derive  scaling  laws  for  the  far-field  period.  Specifically  we  find: 

a.  The  height  of  burst  dependence  is  very  wealQ  for  detonations  in 
the  lower  atmosphere  it  is  the  order  of  jl  ♦  J  where  z  is  the 
height  of  burst  in  kilometers. 

b.  For  yields  small  enough  that  both  weak  shock  theory  and  modified 

Sachs  scaling  are  valid  the  yield  dependence  is  approximately 
1/3 

Y  .  For  very  large  yields  the  situation  is  uncertain. 
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c.  Far  field  periods  ire  found  Co  depend  on  the  altitude  oi  the  tipper 

z  /4II 
o  s 

reflecting  boundary,  zq,  approximately  an  e  '  ,  where  11^  is  tlie 
atmospheric  scale  height.  The  altitude  is  Itself  primarily 
d-tc  rained  bv  t  direction  of  prnpap.it  Ion  relative  to  the  winds 
at  about  50  km. 

3 . 2  I*  ro  pa  gat  ion  Model 

In  this  section  we  derive  the  basic  equation  relating  the 
period  of  the  far-field  high  frequency  acoustic  signal  to  the  mar- 
field  shock  parameters. 

We  assume  that  the  signals  of  interest  in  the  far  field  art 
the  ones  of  largest  amplitude.  Furthermore,  we  assume  that  the  period 
of  these  signals  corresponds  to  the  period  of  the  largest  amplitude 
Fourier  components  of  the  weak  shock  wave.  Clinrncter  Izinp,  the  shock 
wave  as  a  Glasstone  pulse^  with  positive  phase  duration  t+  we  find 
that  the  largest  amplitude  Fourier  component  ius  a  period: 


T  —  2tt  . 


(3-1) 


The  problem  Is  thus  to  calculate  t^,  which  quantity  will 

increase  with  increasing  distance  from  the  detonation.  To  calculate 

the  evolution  of  t,  we  use  Reed's  weak  shock  equations  for  an 

(5) 

inhomogeneous  atmosphere: 
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whore  au<l  P(r)  are  the  values  of  relative  overpressure, 

positive  phase  duration,  and  anblcnt  pressure  at  a  distance  r  from  the 

£P 

dc tom t ion.  ~ — r,  t,  and  P(r  )  are  the  sane  quantities  at  a  smaller 
* tr  J  +  o 

o 

distance  r  ,  y  is  the  ratio  of  specific  heats,  and  c  Is  the  sound  speed. 

These  equations  arc  written  for  a  spherical  shockwave.  In 
the  real  atmosphere  the  shock  wave  is  refracted  and  becomes  non-sphcricnl . 
However,  hydrodynamic  calculations  at  the  Air  Force  Weapons  Laboratory 
have  shown  that  to  a  very  high  degree  of  approximation,  for  an  exponential 
atmosphere  each  portion  of  the  shock  front  can  be  regarded  as  propagating 
independently  with  its  own  radius  of  curvature. The  equations  arc 
also  written  for  an  N  wave  pulse  form.  Groves  has  examined  the  difference 
between  N  wave  and  Glasstonc  pulse  propagation  laws.  He  finds  the 
difference  to  be  very  small. 

The  dependence  of  t+  on  values  at  xQ  may  be  eliminated  by 
taking  the  logarithmic  derivatives  of  Eqs.  (3-2)  and  (3-3)  and  combining 
to  obtain: 


_  Xtir  _£P_ 
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For  an  exponential  atmosphere  this  becomes: 
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where 


n  *  r 
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3r 
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(3-6) 


and  tl  is  Che  effective  scale  height  In  the  direction  of  interest  from 
Che  burst.  The  negative  and  positive  signs  in  hq.  (3-3)  correspond  to 
upward  and  downward  shock  propagation,  respectively. 

Applying  Eq.  (3-5)  at  the  starting  radius  r^  above  the  burst 
and  substituting  back  in  Eq.  (3-2)  yields: 


ty(r) 
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where  Pfa  is  the  ambient  pressure  at  the  burst  altitude  z. 

The  range  dependence  of  t+  is  contained  in  the  Integral: 
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(3-8) 


wtiich  Is  to  be  done  over  a  propagation  path  such  as  illustrated  in 
Fig.  3-1. 

It  is  convenient  to  rewrite  Eq.  (3-8)  as: 


L 


i-l 


where  the  index  £.  corresponds  to  each  straight  line  segment  as  shown  in 
Fig.  3-1,  and  L  is  a  large  number  for  the  far-fiold  disturbance.  We  now 
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turn  Co  Che  Cask  of  sisaning  Che  series  in  Eq.  (3-9).  Assuming  that  the 
radius  of  curvature  is  unaltered  upon  reflection  we  obtain  for  each 

Cera: 
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where  11^  is  the  atmospheric  scale  height  in  tlie  vertical  direction  and 
Che  exponential  integrals  E^  and  arc  defined  ns: 
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El(y)  "  |  t~~  dt  • 
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For  large  y  Eqs.  (3-13)  and  ( 3— 14 )  can  be  approximated  as: 
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On  the  assumption  that  — —  >>  1  we  nay  then  neglect  the  second  terns 

s 

in  Eqs.  (3-11)  and  (3-12)  since  these  are  saaller  than  the  first  teres 

(zo-  z) 

by  a  factor  of  exp - -jjj- - .  Physically  this  is  related  to  the  fact 

8 

that  eost  of  the  contribution  to  the  integral  I,  and  hence  most  of  the 
positive  phase  lengthening,  occurs  at  the  upper  reflecting  boundary 

P^~ 

where  the  quantity  p’(r~  occur  ing  in  Eq.  (3-8)  is  largest. 

Using  the  asvnptotic  fores,  Eqs.  (3-15)  and  (3-16),  for  the 
remaining  exponential  integrals  which  are  functions  of  zq  we  obtain 
the  following  expression  for  1: 
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The  summation  can  be  further  approximated  for  z^  >  z 


(3-18) 
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where  T  is  Euler’s  constant  (.57721)  and  Bp  —  are  Bernoulli 
numbers  (1/6,  1/30,  — ). 

the  first  term  in  this  last  expression  comes  from  I ^ ;  henre 
a  somewhat  better  approximation  than  (3-17)  is: 
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where  m  is  the  number  of  reflections  occurring  at  the  upper  boundary. 

Combining  Kqs.  (3-1),  (3-7),  (3-8),  and  (3-19)  we  finally 

obtain: 
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3.3  Starting  Values  for  Weak  Shock  Parameters 


AP 

It  now  remains  to  determine  the  shock  wave  parameters 

and  n(rQ)  at  rQ.  To  do  this  we  use  modified  Sachs  scaling.  According 
to  modified  Sachs  scaling  the  relative  overpressure,  in  an  inhomogeneous 
atmosphere,  at  a  distance  r  from  an  explosion  where  the  ambient  pressure 
is  P(r),  is  just  the  same  as  if  the  explosion  had  occurred  in  a 
homogeneous  atmosphere  witli  ambient  pressure  P(r).  That  is: 


AP 

P(r) 


(3-21) 
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where  Y  is  the  explosion  yield  and  f  is  the  function  which  gives  the 
relative  overpressure  in  a  homogeneous  atmosphere. 

Modified  Sachs  scaling  appears  to  be  much  better  for  shocks 
which  propagate  downward  in  an  exponential  atmosphere  than  for  those 
which  propagate  upward.  Nevertheless,  modified  Sachs  scaling  is 
superior  to  ordinary  Sachs  scaling  (where  P(r)  on  the  right  side  of 
Eq.  (3-21)  is  replaced  by  P^)  ar.d  it  is  our  hope  Lhat  by  using  it  we 
can  obtain  some  additional  accuracy  in  our  calculations. 

The  accuracy  of  modified  Sachs  scaling  is  best  illustrated 
by  a  particular  example.  In  Fig.  3-2  we  have  compared  the  predic¬ 
tions  of  modified  Sachs  scaling  for  the  relative  overpressure  of  the 
upward  going  shock  with  the  results  of  SAP  and  SHELL  calculations  for 

the  shock  propagation  at  45°  from  the  horizontal  due  to  a  4  Mt  iso- 

(13) 

thermal  sphere  at  sea  level.  The  SHELL  calculation  gives  somewhat 

higher  relative  overpressures  than  does  SAP  because  it  includes  the 
effects  of  fireball  rise  which  SAP  does  not.  The  rising  fireball, 
particularly  for  large  yields,  prevents  the  shock  from  relieving 
backwards.  At  very  small  angles  from  the  horizontal  the  differences 
between  SHELL  and  SAP  are  negligible. 

The  modified  Sachs  scaling  curve  in  Fig.  3-2  has  been 
constructed  under  the  following  assumptions: 

a.  The  4  Mt  isothermal  sphere  is  assumed  to  correspond  to  a  4/. 7  »  5.72 
Mt  detonation.  This  is  based  on  the  value  of  .7  for  the  effective 
yield  for  weak  shock  production  from  nuclear  explosions  in  real 

air  relative  to  ideal  air  given  in  Ref.  14. 

b.  The  1962  Standard  atmosphere  is  used  for  the  ambient  pressure. 

c.  The  homogeneous  atmosphere  calculations  of  Lehto  and  Larson  are 

(14) 

used  to  provide  the  function  f  occurring  in  Eq.  (3-21). 

For  reasons  which  appear  below  we  will  choose  rQ  to  be  the 
AP 

distance  at  which  °y  «  .2.  It  can  be  seen  from  Fig.  3-2  that  for 

o' 
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FIGUR1.  3- 


RELAT1VE  OVER  PRESSURE  ~ 


COMPARISON  OF  SHELL  AND  SAP  CODES  AT  45°  WITH  MODIFIED 
SACHS  SCALING  FOR  4  Mt  SEA  LEVEL  ISOTHERMAL  SPHERE 
(After  Reference  13) 
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this  particular  case  modified  Sachs  scaling  overestimates  rQ  by  perhaps 
5-10%  compared  to  the  SAP  calculation.  The  function  n(rQ)  which  is  the 
reciprocal  of  the  slope  of  the  curve  in  Fig.  3-2  appears  to  be  a  smaller 
negative  number  than  given  by  the  SAP  result.  Quantitative  comparison 
of  the  modified  Sachs  scaling  result  with  the  SHELL  calculation  cannot 
be  done  using  Fig.  3-2  for  directions  of  propagation  other  than  45°. 
However,  for  angles  less  than  45°  the  agreement  between  SHELL  and 
modified  Sachs  scaling  will  be  better  than  shown  in  Fig.  3-2. 

In  any  case  we  would  expect  modified  Sachs  scaling  to  become 


AP 


more  accurate  in  predicting  the  radius  at  which  ^ 

o 

yield  becomes  small.  The  reasons  for  this  are: 


.2  when  the 


1/3 

a.  rQ  a  Y  and  hence  the  atmosphere  encountered  by  the  shock  in 
propagating  to  r  is  more  nearly  homogeneous  for  small  Y. 


b.  The  interaction  of  the  rising  fireball  with  the  upward  going  shock 

(13) 


is  less  important  for  small  yields. 


-r/H 


Returning  to  Eq.  (3-21)  and  putting  P(r)  =  Pb  e  we 
calculate  n(r) 


n(r) 


AP 

P(r) 
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where  f'  is  the  derivative  of  f  with  respect  to  its  argument  and  n^  is 
the  value  of  n  for  a  homogeneous  atmosphere  (H  -*•  ») .  The  quantity  n^ 
is  plotted  vs.  relative  overpressure  in  Fig.  3-3  which  is  based  on 
Ref.  14. 

AP 

We  take  - ^ -y  =  .2  for  which,  according  to  Fig.  3-3, 
o 

n^  ■  -  3/2.  Eq.  (3-22)  then  gives: 
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0-23) 


n(rQ) 


The  radius  rQ  is  given  by  nod  if led  Sachs  scaling  as: 

1/3 

km  ,  (3-24) 

where  Y  is  in  kT  and  P(ro)  is  in  atmospheres.  The  constant  appearing 
in  Eq.  (3-24)  is  obtained  from  Ref.  14. 

AP 

Substitution  of  Eq.  (3-23)  in  Eq.  (3-20)  with  — ■  .2 
gives:  0 
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where  according  to  (3-24): 
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~rt/3H  I  v  r'”  i/3 

rQ  e  -  .53 1  1  -  .53  e  *  Yi/J 
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since  P^  is  measured  in  atmospheres. 

Combining  Eqs.  (3-25)  and  (3-26)  we  obtain  the  basic  equation 
relating  far-field  periods  to  detonation  and  propagation  path  parameters: 
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or  using  the  large  argument  approximation  given  by  tq.  (3-16)  for 
fz  -  5t 
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Equations  (3-28)  and  (3-29)  provide  the  basis  for  the 
discussion  of  scaling  laws  which  follows. 
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3.4  Scaling  Laws  for  the  Far  Field  Period 


3.4.1  Range  Scaling 

The  range  dependence  in  Eqs.  (3-28)  and  (3-29)  is  contained 
in  the  factor  V  1  +  •  Since  n  is  the  number  of  reflections  from 

f  '  ■  ■ 

_  1  R 

the  upper  boundary  this  factor  is  approximately  *  +  2  in  d  *  wbcrc 

K  is  the  range  and  d  is  the  distance  between  ground  reflections.  Based 
on  ray  tracing  calculations  for  infrasonic  propagation  in  a  realistic 
atmosphere^)  a  typical  value  of  d  would  seem  to  be  about  200  km.  The 
~  1  R 

factor  </  1  ♦  2  *'n  d  M0U*d  then  be  equal  to  1.47  at  2,000  km,  2.62  at 

5,000  km,  and  1.72  at  10,000  km;  that  is,  an  Increase  of  17Z  between 

2,000  and  10,000  km. 

Actually,  this  is  the  largest  range  dependence  we  would 
expect  since  processes  we  have  not  considered  would  tend  to  inhibit 
the  increase  of  period  with  range.  These  processes  operate  in  two 
ways:  (a)  They  decrease  the  weak  ahock  amplitude,  for  example  by  pulse 
splitting  due  to  wind  ducts,  and  this  decreases  the  importance  of  non¬ 
linear  effects  of  which  period  lengthening  is  a  manifestation. 

(b)  .Nonlinear  effects  are  also  decreased  in  importance  by  erosion  of 
the  weak  shock  front,  that  is,  by  an  increase  in  the  rise  time  at  the 
front.  This  could  be  caused,  for  example,  by  ambient  turbulence  parti¬ 
cularly  at  higher  altitudes  where  there  is  large  wind  shear.  If  the 

shock  front  is  sufficiently  eroded  the  shock  will  no  longer  satisfy 
the  Rankin e-Hugoniot  relations  on  which  Eqs.  (3-2)  and  (3-3)  are 
based  and  we  would  expect  pulse  lengthening  to  be  less  rapid  than 
these  equations  predict. 

In  summary  the  range  dependence  should  be  very  weak. 
Furthermore,  the  exact  dependence  is  probably  affected  by  meteoro¬ 
logical  conditions  along  the  propagation  path. 
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3.4.2  Height  of  Burst  Scaling 

The  dependence  of  far-field  period  on  height  of  burst  z 
fc”  detonations  Ir  thf  *ower  atmosphere  is  given  by  Eq.  (3-29): 


z/12  li 

s 


1  + 


45.7 


(3-30) 


where  z  is  measured  in  kilometers  and  where  we  have  set  the  scale 
height  II  ■  7  km  and  the  altitude  of  the  upper  boundary  zq  ■  50  km. 

The  origin  of  this  weak  height  of  burst  dependence  can  be 
seen  by  inspection  of  Eq.  (3-29).  The  dependence  on  height  of  burst 

i/3Hs  -1/3 

occurs  in  two  factors:  (a)  e  P,  ,  which  is  due  to  the  fact 

n 


that  the  initial  positive  phase  duration  t°  increases  with  increasing 
height  of  burst  according  to  hydrodynamic  (Sachs)  scaling. 
iz  -z  i 


(b) 


which  represents  the  increase  in 


positive  phase  dura- 


*V  z  -5 

Y  o  o 

at  ion  from  the  initial  value  t+  to  the  value  at  the  upper  reflecting 
boundary.  This  factor  is  a  decreasing  function  of  height  of  burst 
since  the  higher  the  burst  altitude  the  shorter  the  path  length  over 
which  the  positive  phase  duration  can  increase.  As  shown  by  Eq.  (3-30) 
these  two  effects  nearly  cancel  for  detonations  in  the  lower  atmosphere. 


3.4.3  Yield  Dependence 

According  to  Eq.  (3-26)  <  II  <  H  for 
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=  2.3  Mt, 
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for  H  =  7  km  and  P.  *  1. 
s  b 
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For  sea  level  detonations  of  yield  much  less  than  2.3  Mt 
r#  <<  H  and  Eq.  (3-28)  gives: 

T  a  Y1/3  (3-32) 

r  /1 211 

For  larger  yields  the  factor  e  occurring  in  Eq.  (3-28) 

1/3 

would  make  T  increase  more  rapidly  than  Y  .  Before  this  factor  can 
become  effective,  however,  several  other  effects  occur.  First,  the 
upward  going  shock  becomes  strong.  The  maximun  value  of  rQ  allowed 
by  Eq.  (3-26)  is  rQ  ■  3H.  The  reason  for  this  is  that  according  to 
modified  Sachs  scaling  the  upward  going  shock  never  reaches  a  value 
AP 

of  777 — r-  as  small  as  .2  for  r  >  3H.  We  conclude  that  if  weak  shock 
f(r)  o 

o  • 

theory  is  to  be  valid  tq  must  be  considerably  less  than  3H  and  hence 

r  /12H 

the  factor  e  would  not  be  very  important. 

Second,  as  noted  by  Greene  and  Whitaker  in  Ref.  13,  and 
as  demonstrated  by  the  comparison  between  SAP  and  SHELL  calculations 
in  Fig.  3-2,  the  rising  fireball  from  large  yield  detonations  interacts 
with  the  upward  going  shock  and  prevents  it  from  relieving  backwards. 
The  effect  of  this  phenomenon  should  be  to  decrease  the  near-fleld 
positive  phase  duration  and  hence  to  decrease  the  far-field  period. 

ro/12H 

This  would  offset  and  conceivably  could  dominate  the  factor  of  e 

3. A. 4  Dependence  on  High  Altitude  Winds 

When  the  signal  propagation  is  in  the  direction  of  the 
winds  near  50  km  (downwind)  we  expect  the  reflection  altitude  (zq)  to 
be  lower  than  when  the  winds  are  perpendicular  (crosswind)  or  opposite 
(upwind)  the  direction  of  signal  propagation.  The  reflection  altitude 

z  /AH 

occurs  in  Eq.  (3-29)  in  the  factor  e  0  s.  Therefore,  if  Hg  is  taken 
as  7  km,  a  10  km  change  in  z  will  produce  aLnut  a  40%  change  in  T. 
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This  effect  would  appear  to  be  anong  the  most  significant  sources  of 
variation  in  the  far-field  period  for  a  given  yield  and  height  of 
burst.  Note  that  the  far-field  periods  will  be  different  in  different 
directions  for  the  same  detonation. 

In  the  case  of  upwind  propagation,  ray  tracing  calculations 
indicate  that  reflection  mav  not  occu:  at  the  50  km  level.  There  is, 
however,  a  reflecting  boundary  at  about  100  km  altitude.  Applying 

Eq.  (3-29)  for  a  sea  level  detonation  the  period  of  signals  propagating 


inr./4ns 

in  this  channel  is  found  to  be  a  factor  — — - — 

*/Too 


/50 

50/4(1 

s 


4.2  times 


e 

as  large  as  for  signals  propagating  in  the  50  km  channel. 

Finally  we  return  to  Eq.  (3-28)  to  obtain  a  numerical 
estimate  for  the  far-field  periods.  According  to  what  has  been  said 
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regarding  yield  scaling  ve  neglect  the  factors  e  c  .We  take  the 
sound  speed  c  •  .31  km/sec  and  the  scale  height  II  ■  7  km.  We  treat 
the  case  of  a  burst  near  sea  level  (z  a  0).  Results  for  other  burst 
heights  can  be  accomplished  by  means  of  Eq .  (3-30).  Equation  (3-28) 
becomes: 
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where  z  is  in  kilometers. 

o  „ 

T 

In  Fig.  3-4  we  plot  .  vs  z  ,  l'or  z  between  40  and 
r  Y 1  /  3  o’  o 

f:0  km  and  for  m  *  1,2,10,  and  25.  The  value  m=l  neglects  pulse 
lengthening  effects  beyond  the  first  reflection  at  the  upper  boundary. 
The  values  m*2,10,  and  25  would  correspond  to  weak  shock  propagation, 
without  pulse  splitting,  to  larger  distances.  For  a  skip  distance 
d  ■  200  km  the  corresponding  ranges  would  be  400,  2000,  and  5000  km. 

Figure  3-4  predicts  the  far-field  period  to  be  the 
order  of  10  sec  for  a  lkT  low  altitude  detonation.  A  more  precise 
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estimate  would  require  detailed  information  concerning  the  altitude  of 
the  upper  reflecting  boundary  as  well  as  the  approximate  number  of 
transits  between  this  boundary  and  the  ground  for  which  pulse  lenthon- 
ing  effects  are  important.  Determination  of  the  reflection  altitude 
can  be  accomplished  by  ray  tracing.  Determination  of  the  appropriate 
m  value  is  more  complex  since  detailed  meteorological  effects  such  ns 
pulse  splitting  by  wind  ducts  must  be  considered.  However,  t he  fnr- 
field  period  is  insensitive  to  the  exact  value  of  m  and  so  only  a 
rough  estimate  n  •  be  given. 


3.5  Conclusions 

Based  on  the  theoretical  model  outlined  in  the  text  we 
have  reached  the  following  conclusions  concerning  the  short  period 
detonation-produced  infrasound  observed  at  large  distances: 


a.  The  period  for  a  1  kT  detonation  is  the  order  of  10  seconds. 

1/3 

b.  For  small  yields  the  period  scales  as  Y  .  For  large  yields  the 
scaling  is  uncertain.  The  interaction  of  the  rising  fireball 
from  high  yield  detonations  with  the  upward  going  shock  needs  to 
be  investigated  before  the  large  yield  scaling  can  be  determined. 

c.  The  range  dependence  of  tar-field  periods  is  extremely  weak,  as 
an  absolute  maximum  about  a  17%  increase  in  going  from  2,000  to 
10,000  km.  Pulse  splitting  by  wind  ducts  and  the  effects  of 
other  small  scale  meteorological  phenomena  are  expected  to  make  the 
actual  dependence  oven  weaker  than  the  above  value. 


d. 


The  period  dependence  or  height  of  burst 

The  dependence  is  given  approximately  as: 
z  is  the  burst  height  in  km. 


is  al 
T  a 


so  extremely  weak. 

1  +  dbi  wherc 


e.  Probably  the  most  significant  variation  in  far-field  period,  other 
than  that  caused  hy  variations  in  yield,  is  caused  by  the  direction 
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of  signal  propagation  relative  to  the  winds  at  about  SO  km  alti¬ 
tude.  For  downwind  propagation  the  signal  reflection  occurs  at  a 
lower  altitude  than  for  cross  wine,  or  upwind  propagation.  The 

dependence  of  period  on  reflection  altitude  is  approximately  given 
z  /28 

by:  T  a  e  ,  where  z  is  the  reflection  altitude  in  kilometers. 

o 

A  10  km  variation  in  zq  produces  about  a  40%  variation  in  T.  The 
altitude  z^  needs  to  be  determined  by  ray  tracing  for  the  appro¬ 
priate  atmospheric  conditions  and  direction  of  propagation. 
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4.  GENERATION  OF  RAYLEIGH  WAVES  BY 
ATMOSPHERIC  NUCLEAR  DETONATIONS 


4 . 1  Introduction 

This  section  is  an  extension  of  work  contained  in  a  previous 
report^"*^  and  which  will  also  be  presented  elsewhere .  For  convenience 

Ref.  16  has  been  included  as  Appendix  B  of  this  report.  In  this  earlier 
work  contained  in  Appendix  B  we  calculated  the  source  strength  for  Rayleigh 
wave  excitation,  which  is  defined  below,  for  intermediate  altitude  detona¬ 
tions.  These  previous  calculations  were  limited  to  burst  heights  between 
about  30  and  100  km  altitude. 

In  the  present  work  we  extend  these  calculations  to  lower 
burst  altitudes.  Specifically  in  the  case  of  a  10  kT  detonation  wc  r 
have  results  as  low  as  .5  km  altitude  and  for  a  1  Mt  detonation  as  low 
as  2  km  altitude. 

The  extension  to  lower  altitudes  is  not  trivial.  As  discussed 
below,  we  have  had  to  develop  different  methods  of  calculation  than 
were  used  for  detonations  above  30  km  altitude. 

Our  definition  of  the  source  strength  for  Rayleigh  wave  exci¬ 
tation  is  based  on  the  work  of  Toksdz  and  Ben-Menahem.  They  have 

considered  the  far-field  Rayleigh  wave  amplitudes  excited  by  a  point 
source  in  a  homogeneous  atmosphere.  For  three-dimensional  propagation 
in  the  absence  of  nonlinear  effects,  the  source  parameters  contribute 
to  the  amplitudes  almost  entirely  through  a  factor 


|  L(u) j 


oo 


-iwt 

e 


dt 


(4-1) 


which  is  the  Fourier  transform  of  the  point  source  pressure  fluctuation 

6p(t).  In  the  case  of  linear  pressure  pulse  propagation,  5p(t)  a  p^(t)  v i 

where  <5p  (t)  is  the  time  dependent  pressure  pulse  on  the  ground  a  distance 
k 

r  beneath  the  source. 
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Even  in  the  case  where  the  pressure  pulse  propagation  is  non¬ 
linear,  because  of  pulse  lengthening  and  dissipative  effects  we  expect 
that  the  source  strength  for  Rayleigh  wave  excitation  will  be  propor¬ 
tional  to  the  Fourier  transform  of  5p  (t)  r.  The  reason  is  simply  that 

E 

the  source  strength  can  only  depend  on  parameters  at  the  ground  and  not 
on  the  previous  history  of  the  pulse  propagation. 

For  reasons  which  will  become  evident,  we  therefore  define 
the  dimensionless  source  strength  to  be 


S(t) 


2nr 
t  p  li 

g 


CD 

I  c$p (t)  e“iu)t  dt 
J 
0 


(4-2) 


where  t  ■  2tt/u)  is  the  Rayleigh  wave  period  of  interest,  p  is  the  sea 

g 

level  pressure,  and  h  is  the  atmospheric  scale  height.  We  assume  a 
Glasstone  pulse^  at  the  ground 


6p  (t)  =  Ap 

O 


(4-3) 


where  Ap  is  the  peak  overpressure  and  t+  is  the  positive  phase  duration. 

We  then  obtain  the  following  result  from  Eqs.  (4-2)  and  (4-3): 


S(t) 


Ap  r 


1  +  1  TT~ 

2  ' 

p  h 

4 

g 

(4-4) 


Equation  (4-4),  which  is  identical  to  the  definition  used  in 
our  previous  work  (Eq.  (B— 19)  of  Appendix  B),  is  the  basis  for  the  subse¬ 
quent  analysis.  In  deriving  Eq.  (4  —A )  we  have  not  considered  details 
of  Rayleigh  wave  propagation,  so  that  we  are  unable  to  provide  absolute 
values  of  far-field  Rayleigh  wave  amplitudes  (or  phases).  The  present 
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analysis  is  directly  applicable  only  to  estimating  relative  yields  or 
heights  of  burst  for  situations  in  which  source  and  receiver  locations 
are  the  same  for  two  or  more  detonations.  However,  the  shock  parameters 
at  the  earth's  surface,  which  are  calculated  in  the  theory,  are  suitable 
as  inputs  to  existing  Rayleigh  wave  propagation  models  such  as  those 
of  Harkrider  and  Flinn^*^  and  Wickel  and  Whitaker . 

In  Section  4.2  we  calculate  Ap  and  t  as  functions  of  yield 
and  height  of  burst.  The  fact  tliat  the  earth's  atmosphere  is  inhomo¬ 
geneous  plays  a  major  role  in  determining  these  shock  parameters,  at 
least  when  the  height  of  burst  is  comparable  to  a  scale  height.  In 
order  to  treat  this  effect  we  use  the  numerical  shock  propagation  calcu¬ 
lations  of  l.utzky  and  Lehto^^,  for  an  exponential  atmosphere,  in 
conjunction  with  weak  shock  theory.  For  burst  heights  of  much  less 
than  a  scale  height  the  atmosphere  may  be  treated  as  homogeneous  but 

the  shock  becomes  strong,  in  this  case  we  use  the  problem  M  results 

(4) 

to  determine  the  shock  parameters. 

In  Section  4.3  these  values  of  Ap  and  t+  are  used  in 
bq.  (4-4)  to  determine  the  source  strength  for  Rayleigh  wave  excita¬ 
tion  as  a  function  of  yield  and  height  of  burst.  The  dependence  on 
detonation  parameters  is  quite  complex,  but  in  general  we  find: 

a.  The  dependence  on  yield  is  stronger  for  lower  altitude  detona¬ 
tions.  For  example,  for  a  Rayleigh  wave  period  of  20  sec  and 

for  yields  between  500  kT  and  1  Mt  the  source  strength  varies 
.  95 

as  Y‘  when  the  burst  height  is  2  km.  When  the  burst  height 
6  3 

is  10  km  SaY*  for  the  same  Rayleigh  wave  period  and  over 
the  same  yield  range. 

b.  The  variation  of  source  strength  with  height  of  burst  can  be 
quite  significant.  For  example,  for  a  Rayleigh  wave  period  of 
20  sec  and  a  yield  of  100  kT  the  source  strength  increases  by 
a  factor  of  about  2.3  as  the  height  of  burst  is  raised  from 

4  kin  to  20  km.  For  this  yield  and  Rayleigh  wave  period,  S  then 
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continues  ro  increase  slowly  reaching  a  maximum  (about  3  times  the 
value  at  4  km)  at  65  km.  S  then  decreases  with  increasing  height 
of  burst. 

4  . 2  Calculation  of  Shock  Parameters 

We  use  three  different  methods  of  calculating  the  quantities 
Ap  and  t  required  by  bq.  (4-4).  The  method  used  depends  on  the  burst 
height  .is  outlined  in  the  following  discussion. 

(a)  As  in  our  previous  work,  contained  in  Appendix  B,  we  base 
our  analysis  on  the  numerical  calculations  of  Lutzky  and  Lehto^2^  for 
the  case  of  intermediate  altitude  detonations.  These  one-dimensional 
calculations  are  for  downward  shock  propagation  in  an  ideal  gas, 
exponential  atmosphere  due  to  a  point  explosion. 

The  Lutzky  and  Lehto  calculations  are  extrapolated  to 
smaller  overpressures  on  the  ground  than  were  treated  numerically  by 
means  of  weak  shock  theory.  The  basic  equations  for  this  procedure, 
which  are  derived  in  Appendix  B, are  as  follows: 

t+  -  £  8<x)  (4-5) 


^  j  f(x)  e'r/2h  ,  (4-6) 

Pg  r 

where  t+  and  Ap  are  the  positive  phase  duration  and  overpressure  at  the 
earth's  surface  from  a  detonation  at  altitude  r,  h  is  the  atmospheric 
scale  height,  c  is  the  speed  of  sound,  and  p^  is  the  sea  level  pressure. 
The  functions  g(x)  and  f(x)  are  given  by: 


g(x)  =  - 


l±L 

2y 


(4-7) 
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(4-8) 


f(x) 


where  y  is  the  ratio  of  specific  heats, 

at  a  distance  r  <  r  beneath  the  burst, 

o  - 


“P, 


is  the  relative  overpressure 


and  f.j  is  the  exponential  integral 


r  -t 

C 

t.  (x)  *  j  —  dt  .  Tlie  quantities  x  and  n  are  starting  values  obtained 
x 

from  the  numerical  calculations: 


x 


(4-9) 


n 


(4-10) 


Equations  (4-5)  and  (4-7)  are  equivalent  to  Eq.  ( li— 1 0)  of 
Appendix  Bfand  Eqs.  (4-6)  and  (4-8)  are  equivalent  to  Lq.  (15-11)  ol 
Appendix  B. 

The  calculations  of  Lutzky  and  Lehto  are  parameterized 
by  the  quantity: 


M1/3 


hr 


E  I 


(4-11) 


where  p^  is  the  ambient  pressure  at  the  burst  point  and  K  is  the  energy 
released  by  the  point  explosion.  According  to  Ref.  (14)  a  nuclear 
explosion  in  real  air  is  expected  to  be  .7  times  as  effective  in  producing 
a  given  weak  shock  overpressure  at  large  distances  as  the  idealized 
point  source  with  y  =  1.4, 

Eq.  (4-11)  becomes: 


Accordingly,  we  take  =  .7  I.(  and 
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P^(mb) 

h  h  2*1777“ 


h  e 


-r/3h 


P^mb) 
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kf 
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kT 


(4-12) 


values  of 


The  starting  values  x  and  n  obtained  for  the  particular 
used  by  Lutzky  and  Lehto  are  shown  in  Table  4-1. 


TABLE  4-1 


STARTING  VALUES  AT  Ap  /p  -  .1  FOR  APPLICATION 

ro  ro 

OF  THEORY  AS  OBTAINED  FROM  NUMERICAL  CALCULATIONS 


'll 

x 

n 

.05 

3,11 

-4.26 

.1 

2.55 

-3.84 

.2 

1.96 

-3.29 

.5 

1.26 

-2.02 

2.0 

.502 

-1.9 

According  to  Eq.  (4-12)  a  given  yield  Y  and  height  of 
burst  r  determine  o^.  When  the  yield  and  height  of  burst  are  chosen 
to  produce  values  occurring  in  Table  4-1  we  obtain  values  of  x  and  n 
and  hence  determine  t+  and  Ap  by  means  of  Eqs.  (4-4)  -  (4-8).  This  is 
what  we  have  done  for  intermediate  altitude  detonations,  r  >  30  km. 

The  calculational  procedure  is  identical  to  that  of 
Appendix  B  with  two  exceptions:  (i)  We  retain  the  term  £^(r/2h)  in 
Eqs.  (4-7)  and  (4-8)  for  increased  accuracy.  This  term  was  previously 
omitted,  (ii)  The  ambient  pressure  at  the  burst  point,  p^  in  Eq.  (4-11), 
was  previously  determined  by  using  the  1962  Standard  Atmosphere.  For 
the  sake  of  consistency  we  now  determine  p^  by  assuming  that  the 
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atmosphere  is  truly  exponential  as  indicated  in  Eq.  (4-12).  These 

differences  produce  no  significant  changes  in  our  previous  values  of  S 

which  are  given  in  Appendix  B  for  intermediate  altitude  detonations. 

The  above  procedure  does  not  work  at  low  altitudes, 

however.  For  example,  for  a  yield  of  1  Mt  there  are  no  o  values  in 

Table  4-1  corresponding  to  altitudes  between  sea  level  and  30  km.  To 

use  the  method  we  have  outlined  for  low  altitude  detonations,  we  need 

x  and  n  for  additional  o,  values  (particularly  a,  2).  how  these 

h  h 

additional  values  are  obtained  is  described  below. 

(b)  For  detonations  at  low  altitudes  corresponding  to 
>  2  we  extend  Table  4-1  by  adding  an  additional  value  of  o.  ,  namely, 

o  =  a.  corresponding  to  a  homogeneous  atmosphere  (h  =  *)  .  For  this 

"  (14) 

value  of  we  have,  of  course,  x=0  and  from  numerical  calculations 

(see  Fig.  2-3)  we  obtain  n  =  -1.38. 

According  to  Eq.  (4-12)  there  is  no  altitude  or  yield 

(other  than  Y=0)  corresponding  to  =  <*>.  however,  we  add 

to  Table  4-1  in  order  to  interpolate  to  values  between  2  and  <«. 

We  now  describe  how  this  interpolation  is  done. 

In  Fig.  4-1  we  plot  l/o^  +  1  vs.  x  using  the  data  of 

Table  4-1  as  well  as  =  ®,  x  *  0.  Note  that  the  plot,  which  is  semi- 

logarithmic,  is  very  nearly  a  straight  line,  i.e.,  in(l/0j  +  1)  a  x. 

Figure  4-2  shows  a  plot  of  n  vs.  x  using  the  same  x  values,  again  the 

plot  is  very  nearly  a  straight  line,  i.e.,  n  n^  +  ax,  where  nQ  and  a 

are  constants. 

Based  on  the  empirical  evidence  contained  in  Figs.  4-1 
and  4-2  we  use  a  linear  interpolation  of  Rn(l/o^  +  1)  vs.  x  to  obtain 
x  given  o^.  We  then  use  a  linear  interpolation  of  n  vs.  x  to  obtain 
the  corresponding  value  of  n.  This  procedure  works  not  only  for 
between  ®  and  2  but  also  for  example  between  .3  and  2.  Thus,  Eqs.  (4-5) 
through  (4-8)  can  now  be  used  for  arbitrary  yield  and  height  of  burst 
provided  that  the  shock  at  the  ground  is  weak.  When  the  shock  at  the 
ground  is  strong  we  must  use  still  a  different  procedure  as  described 
below. 
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c. .  When  the  shock  at  the  ground  is  strong  (Ap/p  a.  1) 

Kqs.  (4-5)  -  (4-8)  are  no  Longer  valid.  This  situation  arises,  of 

course,  for  low  altitude  detonations.  In  general,  the  height  of  burst 

is  the  order  of,  or  less  than,  an  atmospheric  scale  height.  This 

suggests  that  values  of  the  shock  parameters  obtained  from  numerical 

calculations  with  a  homogeneous  atmosphere  will  be  reasonably  accurate. 

(4) 

We  have  used  Ap  and  t+  values  as  computed  in  problem  M 

with  the  ambient  pressure  at  the  homogeneous  atmosphere  taken  as  p  , 

the  sea  level  pressure.  Scaling  to  sea  level  rather  than  burst  point 

conditions  corresponds  to  modified  Sachs  scaling.  This  scaling  has 

been  shown  to  be  verv  accurate  in  predicting  the  overpressure  Ap. 

The  validity  of  modified  Sachs  scaling  for  strong  shocks  in  predicting 

the  value  of  positive  phase  duration  t+  has  not  been  investigated. 

However,  as  will  be  seen,  this  procedure  does  give  values  of  t+  which 

match  smoothly  to  the  values  given  by  weak  shock  theory  using  the 

cal culational  procedures  (a)  and  (b). 

Figures  4-3  and  4-4  shown  Ap/p^  and  t+  as  functions  of 

r  for  five  different  yields  based  on  this  combination  of  weak  shock 

theory  and  problem  M  results.  We  have  taken  y  »  1.4,  h  *  7  km,  and 

c  -  .34  km/sec,  the  value  used  in  problem  M,  in  constructing  these 

curves.  In  Fig.  4-4  when  neither  weak  shock  theory  nor  the  problem  M 

results  can  strictly  apply  (Af/p  :  1,  r  ;  h)  we  have  used  a  best  fit 

ft 

by  eye.  The  deviations  from  this  best  fit  which  arise  from  applying 
weak  shock  theory  at  too  large  a  value  of  Ap/p  and  applying  the 
problem  M  results  at  too  large  a  value  of  r  are  indicated  by  dashed 
curves  for  the  10  kT  case.  In  all  cases  the  region  fit  by  eye  (or 
"extrapolated")  is  indicated.  For  reasons  described  below  we  do  not 
treat  values  of  Ap/p^  >  1. 

The  solid  curves  in  Figs.  4-3  and  4-4  give  the  values  of 

Ap/p  and  t+  which  are  used  in  Section  4.3  to  predict  the  source 

8 

strength  for  Rayleigh  wave  excitation. 
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WEAK  SHOCK 


PROBLEM  M 


EXTRAPOLATION 


FIGURE  4-4:  DURATION  OF  POSITIVE  PHASE  AT  THE  EARTH’S  SURFACE  VS.  HEIGHT  OF  BURST 


4.3  Results 


To  obtain  the  source  strength  for  Rayleigh  wave  excitation 

we  luve  substituted  the  values  of  Ap/p  and  t,  shown  in  Fics.  4-3  and 

g  + 

4-4  Into  Lq.  (4-4). 

Figures  4-5  through  4-9  show  the  source  strength  as  a 
function  of  altitude  for  yields  of  10,  100,  500,  1000,  and  5000  kT. 

In  each  figure  the  source  strength  is  shown  for  Rayleigh  wave  periods 
of  10,  20,  and  40  seconds.  At  altitudes  above  perhaps  40  km  the  yield 
referred  to  should  be  considered  to  be  an  effective  yield.  This  may 
be  less  than  the  nominal  yield  due  to  radiative  or  other  energv  loss 
mechanisms. 


The  highest  altitude  shown  on  !'ese  plots  corresponds  either 
to  the  smallest  value  of  at  our  disp  al  (o^  =  .05  in  Table  4-1)  or 
to  an  altitude  of  100  Ian.  Our  reason  f^r  stopping  at  100  km  is  that 
above  this  altitude  the  real  atmosphere  cannot  be  characterized  by  a 
constant  scale  height  h  as  assumed  in  the  theory. 

The  lowest  altitude  shown  has  been  arbitrarily  taken  to 

correspond  to  Ap/p  =  1.  Using  the  problem  M  results  we  could  continue 

the  curves  to  still  lower  altitudes.  However,  in  evety  case  we  find 

that  S  rapidly  goes  to  infinity.  This  is  the  result  of  the  fact  that 

3  2 

for  a  strong  shock  Ap  a  1/r  so  that  Apr  a  1/r  -»  n5  as  r  -»  0.  The 
difficulty  seems  to  be  that,  while  the  shock  description  is  correct, 

S  no  longer  represents  the  source  strength  for  Rayleigh  wave  excitation. 
This  is  not  unexpected  since  the  theory  of  Toksoz  and  Ben-Menahem 
which  we  have  used  in  defining  S  is  based  on  a  linear  model  of  wave 
propagation.  A  different  model  of  the  excitation  itself  is  required 
to  go  to  still  lower  altitudes.  Ideally  this  model,  which  would 
include  effects  such  as  cratering,  would  extend  through  the  "trans¬ 
surface"  regime.  That  is,  it  would  match  smoothly  not  only  to  our 
results  but  also  to  results  obtained  for  contained  underground 
explosions. 


50  - 


FIGURE  4-5  DIMENSIONLESS  SOURCE  STRENGTH  AS  A 
FUNCTION  OF  ALTITUDE  FOR  Y  *  10  KT 
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FIGURE  4-6:  DIMENSIONLESS  SOURCE  STRENGTH  AS  A 
FUNCTION  OF  ALTITUDE  FOR  Y  =  100  KT 


"IGURL  4-7:  DI  ML'»S  IONLESS  SOURCE  STRENGTH  AS  A 

FUNG, ION  OF  ALTITUDE  FOR  Y  =  500  KT 
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.DIMENSIONLESS  SOURCE  STRENGTH  AS  A  FUNCTION  OF  ALTITUDE  FOR 
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In  any  case,  the  curves  in  Figs.  4-5  through  4-9  generally 
exhibit  a  maximum  at  intermediate  altitudes  (the  exceptions  are 
Y  >  500  kT  for  x  =  10  sec).  This  behavior  may  be  understood  in  terms 
of  Eq.  (4-4)  in  the  following  way:  At  altitudes  above  the  altitude  of 
the  maximum  27Tt+  >>  x  and  S  a  Apr,  which  is  a  decreasing  function  of  r. 
For  altitudes  somewhat  below  the  altitude  of  the  maximum  2vt+  <  i  and 
S  a  Apr  t  The  positive  phase  duration  t+  decreases  more  rapidly 

with  decreasing  r  than  Apr  increases  and  hence  S  decreases  with  decreas¬ 
ing  r.  The  situation  is  reversed  at  sufficiently  low  altitudes, 

2 

however,  and  Apr  increases  more  rapidly  with  S  than  t_^  decreases. 

Hence  S  may  start  to  increase  with  decreasing  r.  This  behavior  is 
exhibited,  for  example,  by  the  5  Mt  curve  with  x  =  20  sec  and  r  <  10  km. 

For  large  yields  and  short  periods  (i.e.,  Y  >  500  kT , 
x  =  10  sec)  the  condition  2rt+  >  x  holds  at  all  altitudes  and  hence 
S  'u  Apr  which  is  a  monotonic  decreasing  function  of  r. 

Table  4-2  gives  the  burst  altitude  at  which  the  maximum  occurs 
for  various  yields  and  Rayleigh  wave  periods.  As  shown  in  Table  4-2 
the  altitude  of  the  maximum  increases  with  increasing  period  and 
decreases  with  increasing  yield. 


TABLE  4-2 

HEIGHT  OF  BURST  AT  WHICH 
MAXIMUM  SOURCE  STRENGTH  OCCURS 


Yield,  Y (kT)  -+ 


10 

100 

500 

1000 

5000 

n 

M* 

1 

0 

a 

10  1 

36 

17 

— 

— 

— 

H 

! 

v 
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O 

20  1 

81 

65 

55 

45 

39 

40 

89 

71 

62 

57 
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The  increase  in  source  strength  between  the  value  at  low 
altitudes  (say  at  Ap/p  =  1)  and  the  value  at  the  maximum  can  be  sub- 
stantial,  particularly  for  smaller  yields.  For  a  Rayleigh  wave  period 
t  =  no  sec  the  increase  is  a  factor  of  7.6,  for  10  kT,  2.0  for  100  kT, 
1.9  for  500  kT,  1.2  for  1  Mt,  and  1.0  for  5  Mt. 

We  turn  now  to  a  discussion  of  the  yield  dependence  exhibited 
by  S.  Figure  4-10  shows  the  source  strength  as  a  function  of  yield  for 
burst  heights  of  2,  5,  10,  40,  and  60  km. 

It  is  apparent  that  in  general  S  increases  more  rapidly  with 
Y  for  small  heights  of  burst  than  for  large.  This  behavior  is  brought 
our  more  clearly  in  Table  4-3  where  we  show  the  value  of  n  that  results 
from  assuming  a  dependence  S  a  Yn  over  particular  yield  intervals  and 
heights  of  burst. 


TABLE  4-3 

VALUES  OF  n  ASSUMING  S  a  FOR  t  *  20  SEC 


Height  of  Burst  (km) 


2 

5 

10 

40 

60 

10-100 

.87 

.78 

.83 

.48 

.57 

100-500 

.93 

.89 

.68 

.61 

.61 

500-1000 

.95 

.85 

.63 

.63 

.58 

1000-5000 

.65 

.56 

.50 

.50 

This  functional  dependence  of  S  on  Y  can  be  qualitatively 

understood  in  terms  of  Eq.  (4-4).  For  high  altitudes  we  hive  the 

limiting  case  2ut+  >>  i  and  S  a  Apr.  If  modified  Sachs  scaling  is 

approximately  valid,  which  Fig.  4-3  indicates  is  the  case,  and  if  the 

1/3 

shock  is  weak  when  it  strikes  the  earth's  surface  then  Apr  a  Y  , 
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/t\\ 


1/3 

i.e.,  S  u  Y  .  At  low  altitudes  the  limiting  case  is  2?rt+  <<  t  and 

2  1/3 

S  a  (Apr)  t+.  Since  t+  should  scale  approximately  as  Y  we  would 

expect  S  a  Y,  again  provided  the  shock  is  weak  when  it  strikes  the 

earth's  surface. 

4 .4  Conclusions 

Our  detailed  conclusions  are  contained  in  the  previous 
section  in  the  form  of  plots  of  source  strength  for  Rayleigh  wave 
excitation  vs.  yield  and  height  of  burst.  Some  general  conclusions 
based  on  these  plots  are  as  follows: 

a.  The  source  strength  S  increases  less  rapidly  with  yield  Y  the 
higher  the  burst  altitude.  For  example,  for  a  Rayleigh  wave 

period  of  20  sec  and  for  yields  between  100  and  500  kT  we  find 

93  89 

S  a  Y'  for  detonations  at  2  km  altitude,  S  a  Y*  at  5  km, 

S  a  y*^8  at  10  km,  and  S  a  Y*^  at  40  km. 

b.  For  large  yields  and  short  Rayleigh  wave  periods  t  (for  example, 
Y  >  500  kT  for  t  *  10  sec)  the  source  strength  is  a  monotonic 
decreasing  function  of  height  of  burst.  For  larger  yields  or 
longer  Rayleigh  wave  periods  S  is  maximized  for  a  particular 
height  of  burst.  The  altitude  at  which  the  maximum  occurs 
increases  with  increasing  t  and  decreasing  Y.  For  Y  ■  1  Mt  and 
t  ■  20  sec  the  maximum  occurs  at  about  45  km  altitude. 

c.  The  variation  of  S  with  height  of  burst  can  be  appreciable.  For 
example,  for  t  =  20  sec  and  Y  =  100  kT,  S  increases  by  a  factor 
of  2.3  as  the  height  of  burst  is  raised  from  4  km  to  20  km. 
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5.  REPORT  CONCLUSIONS 


Three  different  kinds  of  acoustic  phenomena  related  to 
nuclear  test  detection  and  diagnostics  have  been  discussed  in  this 
report.  They  are:  (1)  Long  period  (3-10  minutes)  acoustic-gravity 
waves  which  propagate  in  the  real  atmosphere's  counterpart  to  Lamb's 
edge  mode.  (2)  High  frequency  (1-60  second  period)  acoustic  pulses 
which  result  from  the  degeneration  of  the  detonation  blast  wave  into 
a  weak  shock.  (3)  Rayleigh  or  surface  seismic  waves  which  result 
when  the  detonation  shock  wave  strikes  the  earth.  Conclusions  reached 
for  each  of  these  phenomena  are  given  below. 

5.1  Lamb  Mode  Analysis 

(a)  The  shock  front  extending  from  the  ground  to  many 
scale  heights  above  the  ground  is  of  almost  equal  importance  in 
determining  the  far-fie  d  Lamb  mode  amplitude,  because  of  this,  the 
far-field  amplitude  should  be  relatively  independent  of  phenomena 
peculiar  to  any  small  portion  of  the  shock  front.  In  particular, 
ground  effects  such  as  the  presence  of  a  precursor  should  not  play  a 
very  important  role. 

(b)  For  large  yield,  low  altitude  detonations  we  have  found 

a  far-field  amplitude  dependence  which  varies  as  Y 

where  typically  Y*  is  the  order  of  10  Mt  for  a  sea  level  detonation. 

This  dependence  arises  from  a  detailed  consideration  of  the  mechanism 
by  which  the  low  frequency  components  of  the  near-field  (shock)  distur¬ 
bance  become  spatially  separated  from  the  high  frequency  components. 

The  precise  value  of  the  critical  yield  Y*  where  yield-amplitude  pro¬ 
portionality  breaks  down  is  determined  by  atmospheric  conditions  near 
the  burst  as  well  as  the  direction  of  propagation  and  must  be  determined 
by  detailed  numerical  calculations.  Nevertheless,  in  analyzing  the  test 
data  we  have  seen  some  evidence  that  very  large  vields  do  on  occasion 
produce  anomalously  large  far-field  amplitudes. 
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(c)  A  number  of  potentially  important  source  effects  can  be 
explored  within  the  framework  of  the  Lamb  mode  theory.  The  most  impor¬ 
tant  of  these  effects,  which  we  have  not  investigated,  are  probably 
related  to  the  presence  of  a  secondary  shock  front  due  to  ground 
reflection.  We  expect  that  at  low  altitudes  the  effect  of  the  secondary 
shock  in  determining  the  far-field  amplitudes  is  approximately  equal 
to  that  of  the  direct  shock,  but  for  higher  altitude  detonations  inter¬ 
ference  or  destructive  effects  between  the  two  shock  may  occur. 

5. 2  High  Frequency  Acoustic  Periods 

(a)  At  large  distances  the  dominant  acoustic  pulse  period 
for  a  1  kT  detonation  is  che  order  of  10  seconds. 

1/3 

(b)  For  small  yields  the  period  scales  as  (yield)  .  For 
large  yields  the  scaling  is  uncertain.  The  interaction  of  the  rising 
fireball  from  high-yield  detonations  with  the  upward  going  shock  needs 
to  be  investigated  before  the  large  yield  scaling  can  be  determined. 

(c)  The  range  dependence  of  far-field  periods  is  extremely 
weak,  as  an  absolute  maximum  about  a  17%  increase  in  going  from  2,000 
to  10,000  km.  Pulse  splitting  by  wind  ducts  and  the  effects  of  other 
small  scale  meteorological  phenomena  are  expected  to  make  the  actual 
dependence  even  weaker  than  the  above  value. 

(d)  The  period  dependence  on  height  of  burst  is  also 
extremely  weak.  The  dependence  is  given  approximately  as:  T  c* 
where  T  is  the  period  and  s  is  the  burst  height  in  km. 

(e)  Probably  the  most  significant  variation  in  far-field 

period,  other  than  that  caused  by  variations  in  yield,  is  caused  by 

the  direction  of  signal  propagation  relative  to  the  winds  at  about  50  km 

altitude.  For  downwind  propagation  the  signal  reflection  occurs  at  a 

lower  altitude  than  for  crosswind  or  upwind  propagation.  The  dependence 
of  period  on  reflection  altitude  is  approximately  given  by: 
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z  /  28 

T  a  e  ,  where  z  is  the  reflection  altitude  in  kilometers.  A  10  km 

o 

variation  in  z  produces  about  a  40%  variation  in  T.  The  altitude  ?. 

o  o 

needs  to  be  determined  by  ray  tracing  for  the  appropriate  atmospheric 
conditions  and  direction  of  propagation. 

5.3  Rayleigh  Wave  Amplitudes 

Our  detailed  conclusions  are  contained  in  the  text  in  the 
form  of  plots  of  source  strength  for  Rayleigh  wive  excitation  vs.  yield 
and  height  of  burst.  Some  general  conclusions  are  as  follows: 

(a)  The  source  strength  S  increases  less  rapidly  with  vield  Y 
the  higher  the  burst  altitude.  For  example,  for  a  Rayleigh  wave  period 

O'! 

of  20  sec  and  for  yields  between  100  and  500  kT  we  find  S  a  Y’  for 
detonations  at  2  km  altitude,  S  a  Y‘^  at  5  km,  S  a  Y*^  at  10  km, 
and  S  a  Y  at  40  km. 

(b)  For  large  yields  and  short  Rayleigh  wave  periods  t  (for 
example,  Y  >  500  kT  for  x  =  10  sec)  the  source  strength  is  a  monotonic 
decreasing  function  of  height  of  burst.  For  larger  yields  or  longer 
Rayleigh  wave  periods  S  is  maximized  for  a  particular  height  of  burst. 
The  altitude  at  which  the  maximum  occurs  increases  with  increasing  x 
and  decreasing  Y.  For  Y  =  1  Mt  and  x  »  20  sec  the  maximum  occurs  at 
about  45  km  altitude. 

(c)  The  variation  of  S  with  height  of  burst  can  be  appre¬ 
ciable.  For  example,  for  x  ■  20  sec  and  Y  =  100  kT,  S  increases  by  a 
factor  of  2.3  as  the  height  of  burst  is  raised  from  4  km  to  20  km. 

(d)  The  present  analysis  of  the  Rayleigh  wave  excitation 
mechanism  appears  to  break  down  when  the  shock  striking  the  ground  is 
strong.  This  occurs  for  1  kT  at  burst  heights  below  about  .5  km  and 
for  1  Mt  at  burst  heights  below  about  2  km.  A  new  model  needs  to  hi 
formulated  between  the  altitude  where  our  model  no  longer  holds  and 
the  burial  depth  of  a  contained  underground  detonation. 
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AiTiifHMX  A:  l)i:RlVATIO.<  OF  Till:  BASIC  HQUATlONS  FUR  P  Kill)  1CT  IPX 
OF  FAR-F I  ELI)  LAMB  MOI1K  AMIM.ITUiJKS 


A. 1  introduction 

In  this  appendix  we  develop  the  theoretical  basis  for  the 
equations  discussed  In  the  text.  The  analysis  is  heavily  based  on  the 
work  of  Fierce  and  Posey  described  in  Ref.  1,  and,  where  possible,  vi 
use  the  notation  of  this  reference.  However,  because  we  are  primarily 
concerned  with  the  source  modi i  rather  than  with  tin-  proposal  Ion  theory 
we  make  several  simplifying  assumptions.  Specifically: 

(a)  We  neglect  the  effects  of  acoustic  ray  refraction  and  nssimn 
that  all  rays  follow  groat -circle  pat hr  between  source  and 
observation  point.  This  is  not  necessarily  a  good  assumption. 
However,  the  corrections  necessary  to  incorporate  a  realistic 
atmospheric  propagation  model  do  not  affect  the  way  in  which  tin- 
source  model  is  Incorporated  into  tin  tbeorv. 

(b)  On  the  assumption  that  the  various  "sound  speeds"  occurring  In 
the  P  lerco-l’osc-y  theory  are  of  approximately  equal  nvignitud< 

we  do  not  distinguish  them,  in  the  notation  of  Ref.  1  this  Means 
that  the  altitude  dependent  sound  speed  e(z),  the  sound  speed 
averaged  in  altitude  over  the  Lamb  mode  weighting  function  c^  , 
and  the  effective  sound  speed  c^  which  incorporates  effects  such 
as  t  ho  sc  due  to  winds,  are  ail  denoted  by  c.  At  tin  same  time, 
and  this  is  consistent  with  the  above  assumption,  the  major  eft  eel 
of  differences  between  these  quantities,  namely  a  frequency 
dependent  group  velocity,  is  taken  Into  account. 

A . 2  Isolation  of  the  i.anb  Mode  In  the  Xcar-Fleld  Disturbance 

According  to  Pierce  and  I’osey,  it  the  upper  atmosphere  Is 

considered  to  be  bounded  at  some  large  altitude,  z  .  by  a  rigid  c«.  lilt 

n 

then  an  acoustic  disturbance  mnv  be  expanded  as  a  sum  of  guided  mod**;. 


one  of  which  is  the  analogue  of  Lamb's  mode.  That  is: 


AF(r,t)  -  l‘,  (r.t)  +  V  ( r,t ),  (A-l) 

win  rc  At’  is  the  overpressure  of  the  disturbance.  is  the  contribu- 
t Ion  of  the  Lamb  mode,  and  T  represents  the  contribution  of  other 
modes.  The  Lamb  mode  is  characterized  by  an  altitude  dependence  such 
that  I*  can  be  factored  as: 


'•L(r.O 


.  I' 

s 


c-z/yll 


F(r,t), 


(A-2) 


when  I’  is  the  sea-level  pressure,  z  is  the  altitude,  y  is  the  ratio 
of  specific  heats,  and  II  is  the  scale  height.  F(r,t)  in  a  function 
only  of  time  and  of  the  great-circle  distance  r. 

The  height  profiles  of  the  various  modes  contained  in  * 
satisfy  an  orthogonality  relation  such  that: 


U 


(A-3) 


Turning  now  to  the  case  in  which  AP(r,t)  represents  the 

detonation  shock  wave,  we  define  tq  to  be  the  appropriate  horizontal 

distance  for  matching  the  shock  to  the  linear  propagation  theory.  How 

this  distance  is  to  be  determined  is  left  for  subsequent  discussion. 

Suppose  that  at  distance  r  and  time  t  the  shock  extends  between 

o  o 

altitudes  Zj  and  Zy  if  will  be  shown  below  that  for  detonations  in 
the  lower  atmosphere  z^  =  0  (i.c.,  sea  level)  and  z^  »  H  is,  in  fact, 
the  situation  of  interest.  Under  these  circumstances  we  assume  that 
i.q.  (A-3)  is  approximately  correct  with  the  limits,  0,z  replaced  by 

R 
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y-1  z 
Y  11 

Multiplying  Kq.  (A-l)  by  e  ,  integrating  from  z=z^  to 


z*Z2,  and  using  Eq.  (A-2)  we  then  obtain: 


i 

/ 


y-1  z 

AP(r  ,z,t  )  e  y  **  dz 
o  o 


F(r  ,t  )  P 
o  o  s 


/ 


z?  1 ill 


Y  H  , 

e  dz 


y-2  z 

F(r  ,t  )P  -hr  e  1 
o  o  s  y~2 


=  F(r  ,t  )V  , 

o  o  .s  2-y 


(A-4) 


2-y  1  2-y  2 

where  we  have  assumed  — L  rr-  <<  1  and  -  x  rr-  >>  1. 

Y  h  y  II 

AP(ro,z,to)  . 

Putting  - jj -  e  ■  :r(ro,z,to)  where  n  is  the  relative 

s 

overpressure,  we  may  rewrice  Eq.  (A-4) 

rZ2  -  — 

F(to>to)  •  fii1 1  a<ro'z'to)  *  «'*•  <A'5) 

•i/ 

Z1 

We  retain  the  limits  z^,  z ^  rather  than  setting  them  to  0,™  as  a 
reminder  of  the  assumptions  made. 

In  the  Pierce-Posey  theory  the  geometrical  acoustics  approxi¬ 
mation  is  made  that  F  may  be  factored  into  amplitude  and  waveform 
functions: 


F(r  » t  )  -  A(r  )  v  (r  ,t  ) 
o  o  o  o  o 


(A-G) 
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and  similarly  at  some  r, t  corresponding  to  a  larger  distance  and  a 
later  time: 


F(r,t)  -  A(r)  •;  (r,t). 


(A— 7) 


For  Eq.  (A— 7 )  to  be  valid  the  spatial  variation  of  A(r)  must  be  on  a 
scale  large  empared  to  the  disturbance  wavelength.  This  will  be  the 
case  for  the  effect  of  earth  curvature  which  produces  a  dependence: 


A(r) 


(A-8) 


where  rfi  Is  the  earth's  radius  and  where  r  is  again  the  great  circle 
distance. 


A. 3  Propagation  of  the  Lamb  Mode  to  the  Far  Field 

The  waveform  function  £>  is  representable  as  a  Fourier 

Integral: 


v(r,t)  -  2Re  j  w(k)  e’1<UJt'kr)  dk,  (A-9) 

o 

The  quantities  u  and  k  are  related  through  the  Lamb  mode  dispersion 
equation: 


.  •  <=  k-k3  V  W-10) 

where  c  is  taken  as  an  average  sound  speed  and  the  quantity  h^  is 
not  a  function  of  k  and  is  due  to  deviations  of  sound  speed  and  wind 
velocity  from  vertically  averaged  values.  The  second  term  on  the 
right  side  of  Fq.  (A-10)  is  much  smaller  than  the  first. 
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Combining  Eqs.  (A-9)  and  (A-10)  we  obtain: 


at 


+  c 


9r 


+  h 


kk 


0 


(A-U) 


The  last  and  smallest  term  in  Eq.  (A-ll),  h^  — ^ 

3r 

dispersive  effects,  may  be  rewritten  by  replacing 


,  which  represents 
to  lowest  order 


i 

c  at  ' 


Thus: 


I  it  +  ii  +  V  ai.  .  0 

c  Jt  3r  4  .3  u 

c  at 


( A— 1 2) 


Following  Pierce  and  Posey  we  now  introduce  the  parameters: 


t 


r 


r 


(A-13) 

(A-14) 


.0  that  |i.  &  and  |i  .  ^  3l  +  -tt  i±  . 
3t  -r  3r  c  -  4 

at  at  c  ar 


Thus  Kq.  (A-12)  becomes: 


it.  A.  o 

-  -3 

as  at 


( A— 1 5) 


Equation  (A-15)  may  be  solved  to  obtain  ij  at  r,t  (r,t)  from 

<J»(t  )  at  r  [*(t  )  at  r  1  . 
o  o  l  o  o  J 
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(r,t) 


4-  (3l :)in  /  *l  Ohin\^ 

■■CP  •• 


1...  - 


r  ,t  )  dt 
o  o  o 


(A-16) 


where  is  the  Airy  function. 

We  now  define  the  waveform  arrival  time: 


t  *»  t  -  t 
a 


(A-17) 


and  the  characteristic  dispersion  time: 


<d  ■  (3?)1/3 


(A-18) 


in  terms  of  which  Eq.  (A-16)  mr-v  be  written 


00  r 

"  7  f  \  :j-Tr— °  *<to’ro,<ico  <A-J 

/it  t,  ^  d 

d  — J 


In  writing  Eq.  (A-19)  we  have  put  to  ■  tQ  since  tQ  is  being  measured 
a1.  rQ,  the  distance  at  which  the  Lamb  mode  theory  is  presumed  to  just 
► Jgin  to  be  applicable. 

The  far-field  overpressure  at  sea  level  is  given  by  combining 
Eqs.  (A-2),  (A-7),  and  (A-19) 


P. (r, t)  -  Pe  F(r,t) 

L  S 


P  A(r)  *(rft) 
s 


*  P  A( 
s 


r)  f  A.I^- 

^  td  L  ‘l  t 


T  -  t  +  t 


*(tn ,r  )dt  (A-20) 

0  0  o 
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Henceforth  we  will  measure  the  time  t  relative  to  the  arrival 


F(ro,t0) 

time  t  .  Now  ip(r  ,t  )  -  —r~. — r —  ,  where  F(r  ,t  )  is  given  by  Eq.  (A-5)  . 
a  o  o  Air  ,)  o  o 

o 

Thus  we  obtain: 


A(r) 
s  A(ro) 


(2-y)  1 
YH/tT  Td 


(A-21) 


A. 4  The  Shock  Model 

To  proceed  it  is  necessary  to  make  some  additional  assump¬ 
tions  about  the  form  of  it.  We  assume  that  the  detonation  occurs  at 
altitude  z  and  that  the  shock  front  is  spherical  about  the  burst.  For 
this  assumption  to  be  correct  it  is  necessary  that  while  the  shock  is 
strong  and  supersonic  it  encounters  a  nearly  homogeneous  atmosphere. 

When  the  shock  is  weak  and  sonic  the  shock  parameters  may  vary 
drastically  along  the  front,  although  the  front  will  still  be 
spherical. 

The  assumption  of  a  single  spherical  shock  front  is  admittedly 
an  oversimplification.  For  detonations  in  the  lower  atmosphere  there 
will  be  a  second  shock,  due  to  reflection  from  the  ground,  which  will 
follow  the  main  shock.  Near  the  ground  the  primary  shock  and  the 
ground  reflected  shock  may  fuse  and  form  a  Mach  stem.  Very  crudely 
the  reflected  shock  may  be  considered  as  arising  from  a  mirror  detona¬ 
tion  at  -z.  The  time  lag  between  primary  and  reflected  shocks  will  be 
2z 

less  than  -  —  and  the  difference  in  radius  of  curvature  between  the 
c 

two  shocks  will  be  less  than  2z.  For  small  z,  then,  the  reflected 
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shock  follows  immediately  behind  the  prinary  shock  and  has  approxi¬ 
mately  the  same  curvature.  The  time  lag  between  the  two  shocks  will 
be  less  than  a  typical  far-field  quarter  period  (-100  sec)  for  z  <  15  km. 
On  this  basis  we  expect  that  for  low  altitude  detonations  the  effect 
of  the  ground  reflected  shock  is  to  approximately  double  the  far-field 
overpressure.  Thus  the  effective  yield  of  a  low  altitude  detonation 
would  be  essentially  the  same  as  for  a  contect  (z-0)  surface  burst, 
although  this  is  not  exact  since  the  contact  detonation  expends  some 
energy  in  digging  a  crater,  etc.  This  model  must  break  down,  however, 
for  large  heights  of  burst  where  the  time  lag  between  the  two  shocks 
becomes  comparable  with  a  far-field  quarter  period  or  where  the  radius 
of  curvature  of  the  two  shocks  at  the  matching  radius  rQ  becomes  very 
different . 


The  techniques  we  develop  below  should  be  adequate  to  treat 

this  more  complex  shock  structure  in  detail  and  we  hope  to  do  this  in 

the  future.  However,  in  the  present  analysis  we  consider  only  a 

single  spherical  shock  front  and  take  the  ground  reflected  shock  into 

account  by  simply  doubling  the  amplitude  in  Eq.  (A-21) .  The  geometry 

envisioned  is  shown  in  Fig.  A-l. 

Let  the  time  t  be  measured  from  the  arrival  of  the  hori- 
o 

zontal  portion  of  the  shock  front  at  r  .  At  time  t  >  0  the 

o  o 

spherical  shock  radius  is  rQ  +  ctQ,  where  c  is  the  sound  speed,  since 

the  shock  is  assumed  weak  and  near  sonic  at  r  . 

o 

For  the  situation  illustrated  in  Fig.  A-l  we  must  have 


z  <  z2  =  z 


/2  2 

c  t  +  2r  ct  and  z  >  z,  •  larger  of  zero  or 
o  o  o  1 


-  '22 

z  -  let  +  2r  ct  .  The  integration  in  Eq.  (A-21)  occurs  over  the 
y  o  o  o 

region  of  space  and  time  within  which  the  shock  is  located,  that  is 
0  <  tQ  <  ”,  <  z  <  z2>  This  is  equivalent  to  an  integration  over 


0  <  z  < 


(z-i)2  +  r2  -  r 
\  _ o  c 


<  t  <  ».  Thus,  switching  the  order  of 
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SHOCK  GEOMETRY  SEAR  THE  HIRST 


Integral  ion  we  obtain: 


.  p  m  ML  f  d,  r*/>H  _  dt 

*  A(r«»)  i H*  'd  j  'A  2.  2  c 

o  Y<*-*  >  *  *,  -  ro 


*  t  -tl 

A  — —  -(r  .z.t  ), 

l.  Td  i  "  u 


<A-22) 


whore  a  factor  of  2  lias  been  added  to  represent  the  contribution  of 
the  ground  roflec ted  shock. 

We  define  the  local  time  t: 


r  +  ct  -  "  J  Li*  ♦  r^ 

t  -  --  —  - - L_ - «  (A-23) 

c 

where  As  •  (z-z).  The  quantity  t  is  the  tine  elapsed  at  ro»«  since 
the  passage  of  the  sluick  front. 

Now  we  further  assume  that  the  relative  overpressure  at 
has  the  form  of  a  (.lasstonc  pulse  when  written  in  terms  of  the 
local  time: 


!(r  .z. t  )  -  *(r  ,z)  e 

O  O  0 


■t'tt(v,)i 


1  - 


t+(r0’l) 


(A-24) 


for  t  >  0.  Equation  (A-22)  becomes: 
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.  p  Mil  j_  f'f\  1*.'/U**  t-T* 

L  s  A(r  )  j—  llx,  j  t  x,  ci. 

o  y*z  d  J  J  d  d  , 

o  o  — 


n(r  ,z)  c  +  (1  -  x/t.)  dx  dz 

o  + 


(A-25) 


As  shown  by  Fierce  and  Foscy  the  characteristic  dispersion 

time  T'j  is  the  order  of  the  far-flcld  waveform  period  which  is  much 

larger  than  the  near  field  positive  phase  duration  t  .  For  x  comparable 
-x/t+ 

to  x^  the  factor  e  occurring  in  the  integrand  of  Eq.  (A-25)  is 

very  small.  Since  the  dominant  contribution  to  the  x  integration  comes 
from  x  <<  Xj  wc  may  expand  A^  in  a  power  scries  in  —  : 


A 


i 


Putting  x  »  wc  do  the  x  integration,  noting  that : 
N- 


(A-26) 


/ 


o 


(1-x)  dx  -  0, 


(A— 2  7 ) 


J  x  c  X  (1-x)  dx  ■  -1  . 
o 

The  result  is: 


(A-28) 
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A(r)  2-y  1 

*  A(ro>  Htd2 
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/U‘2  *  ro 
ctd 


c-2/yH 
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(r  ,z)  t 
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2 
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(rQ,z)  dz 


(A-29) 


Finally,  using  Lq.  (A-8)  under  the  assuaption  that  r  «  r 

o  c 

we  obtain: 


P 


L 


212lll  _L_ 

\S*  Hi  * 

a 


I 

* 

o 


+ 


c,d 


c-z/yh 


*)  t. 


(rQ,z)  dz 


(A- 30) 


A. 5  Form  of  the  Shock  Parameters  it,  t+ 


We  wish  to  do  the  z  integration  in  Eq.  (A-30)  using  a  simple 
shock  model  for  *t (r  , z)  and  r  (r  , z).  The  shock  aodel  assuaes  that  at 

o  TO 

some  snail  radius  Rq  <  r  about  the  burst  the  shock  parameters  it(RQ), 
t+(R  )  arc  constant  along  the  spherical  shock  front  and  that  the  shock 
is  weak  at  R^  ( - ( )  «1).  The  geoactry  is  shown  in  Fig.  A-2.  It  is 
only  necessary  that  this  assuaption  apply  to  the  portion  of  the  shock 
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FI  CTRL  A-2:  SHOCK  GEOMETRY  AT  SMALL  RADIUS  Rc  WHERE  ATMOSPHERIC 
INHOMOGENLITY  IS  NEGLECTED  AND  AT  A  LARGER  RADIUS 
WHERE  IT  IS  NOT 
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front  that  is  effective  in  exciting  the  Lamb  mode;  that  Is,  it  need 
not  apply,  for  example,  to  the  vertical./  upward  going  portion  of  the 
shock  front. 

Beyond  this  sn.il  1  radius  R  the  shock  Is  a  sawed  to  evolve 

0  (2) 

according  to  the  weak  shock  equations.  In  particular  at  (r  ,z): 


i(r  ,z) 
o 


R  C4.(R«>  /  P(Rn) 

\  _ 2 _ 1_9 _  ,  /  ° 

°  ,  rrm  t+<r.,*)  *  p(vi) 


(*-31) 


r  +  Az 


C+(ro‘z) 


„r  ,z  _  ^  4 

>+i  "<Ra>R  ®  /P(R  )  dR 

t+(Ro)  ^  +  2y  ct+(«o)  / “RrT  IT 


(A-32) 


where  l’(K  )  is  the  ambient  pressure  at  R  ,  P(r  ,  z)  is  the  ambient 
o  o  o 

pressure  at  the  point  r^.z  where  matching  to  the  Lamb  mode  propagation 
model  occurs,  and  i’(R)  Is  the  ambient  pressure  at  a  spherical  radius  R, 

V2  2 

rQ  +  A z  .  Consistent  with  the  assumption  that  Rq  Is  a 

snail  radius  (that  is,  that  It  is  small  compared  with  the  effective 

scale  height  in  the  direction  In  question)  we  take  P(R  )  •  P.  ,  the 

o  b 

-Az 

u 

ambient  pressure  at  the  burst  point.  Then  P(ro»z)  ■  P^  e  and 
Kq.  (A— 31 )  becomes: 


-<ro.Z) 


!T 

O 


2  .  ,  2 
r  +  Az 
o 


t+(vt> 


(*-33) 


where  (except  in  the  case  of  rQ)  we  use  zero  subscripts  to  denote 

values  at  K  . 

o 
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The  li.tegr.il  occurring  in  hq.  (A-32)  may  be  approximated  as 


follows: 


r'°'\ 


"<V  dK 

P(R)  K 


“  /r2  ♦  Az2 

r 


R  CXP 


(K-H) 

o 


■>  •> 
*  ,  * 
r  +  Az 
.  o 


ft ’  "  2 
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/  o 


-  2  ?. 
V  ro  +  ** 
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I  ' - =^=  ir 

n-o  211  r  +  Az 
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R"”1  dK 


(A-3A) 


■»  2  2  1  -  2  2 
\  ro  +  Az  /  En  Vro  +  LZ 


f  M"  -L 

L.  [211/  nn! 


V'r2  *  l«2l  j  f 

V  °  J  j 


Consistent  with  the  assumption  tli.it  Rq  is  small,  i.e. 


r 

R„  <  y  r2  +  Az2  ,  R  <<  211 
o  V  o  o 


V7 


,  S.q.  (A-34)  reduces  to: 
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(A-35) 


Thus  Eq.  (A-32)  becomes: 
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-o'* 


t+(ro.z)  -  t°  1  ♦  6  in 

L  - 
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2.2 
r  +  Az 
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♦  z, 

n»l 


1  i  A* 

nn!  j  211 1 


J 


/  ,  (A-36) 


.  |  n  r\ 

where  we  have  defined  8  •  The  quantity  8  is  nearly  constant 

for  a  weak  slwck.  Since  we  want  Rq  as  small  as  possible  consistent 
with  the  shock  being  weak  wc  take  7q  ■  .1.  Than  6  ■  .67  *  2/3  accord- 

(31 

ing  to  the  problem  N  results  with  y  ■  1.4  and  c  ■  .34  km/sec.  The 

other  parameters  used  in  computing  8  correspond  to  a  distance  of  .9  km 

from  n  1  k'f  explosion  where  the  positive  phase  duration  is  about  1/3  sec. 

Now  in  our  basic  equation,  Eq.  (A-30),  the  shock  parameters 
.  2 

occur  in  the  form  n(r  ,z)  t.(r  ,  z).  Using  Eqs.  (A-33)  and  (A-36)  we 

O  ■TO 

write  this  quantity  in  terns  of  the  value  at  rQ,  z-z,  that  it  the  value 
for  the  horizontally  propagating  portion  of  the  shock: 


’ (r  ,z)  t  (r  ,z) 
o  +  o 


fl(r  •  z)  t  (r  ,z) 
o  ▼  o 


»s 


(A-37) 
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Thus  tq.  (A-30)  becomes 
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(A-38) 


.2  2 
Az  +  r  -  r 

We  now  introduce  the  variable  y  ■  — 1 - 2 - -  and 

J  ci. 


|  Az| 


.  The  y  integration 


0  <  y  <  ®  (Az  >  0)  and 

-V  [-2  A  2 

\/  z  +  r  -  r 
V  o  o 


0  <  v  < 


ci 


(Az  <  0). 
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Equation  (A-38)  becomes: 


where  in  terms  of  the  variable  Az(y): 


We  will  demonstrate  below  that  G(y)  is  a  weak  function  of  Az  and 
hence  of  y  over  the  integration  region  of  interest.  Accordingly,  as  a 
first  approximation,  we  take  G(y)  to  be  given  by  its  value  at  y*0  (Az*0) 
which  is  unity.  The  remaining  integral  for  Az  >  0  is  twice  the  function 
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da  discussed  by  Pierce  and  Posey. 


PP 


|i- )-  f  a;  [-  —  +  a2 1 

I’d!  J  1  l  Td  J 


(2) 


They  have  shown  by  numerical  computation  that  the  first  maximum  of 

PP  [— |  occurs  at  —  =  1.  The  remaining  integral  for  Az  <  0  will  be 
'  \  Td  /  Td 

shown  to  be  unimportant.  Therefore,  for  consideration  of  the  first 
maximum  of  the  far  field  amplitude  we  may  restrict  our  attention  to 


times  t  ~  t,. 

d 


AI  .-1  +  yj 

The  integrand -  ,  which  results  from  setting 

✓  y 

t  =  ,  is  plotted  in  Fig.  A-3.  Maximum  contribution  to  the  integral 

is  seen  to  come  from  the  region  0  <  y  $  2  corresponding  to 


o  <  i  az  i  <  2crd  yr+r 7  CTj.  Since  is  a  quarter  period  observed 

in  the  far  field  a  typical  value  is  about  100  seconds.  For  c  -  .3  km/sec 
we  thus  find  that  | Az |  values  up  to  at  least  60  km  contribute  to  the 
integral  in  Eq.  (A-39) .  This  justifies  the  initial  assumption  made  in 
extracting  the  Lamb  mode  from  the  near-field  pressure  disturbance; 
namely,  that  the  z  integration  extends  to  values  much 


larger  than  a  scale  height,  z 2  >>  h. 

It  remains  for  us  to  justify  the  assumption  made  above  that 
G(y)  -  1  over  the  range  of  Az  of  interest.  First,  however,  we  note 
the  physical  origin  of  the  various  factors  composing  G.  These  factors 
are: 

-Az 

yH 

a.  e  is  the  Lamb  mod(  weighting  factor. 

Az 

2h 

b.  e  is  due  to  the  increase  in  relative  overpressure  tt  with 
altitude. 
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1  + 
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_  ■  ...  Az^  1  /«lz 

v  7*  ' 

/  o  n«l 
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>  *5 


1  +  B  In  ^ 
o 


'J 


is  due  to  the  change 


d. 


in  positive  phase  duration  t+  with  altitude 

/2 


/ 


is  a  geometrical  factor  arising  from  Lhc 


V,  + V1*?1 


assumption  of  a  spherical  shock  front. 


The  relatively  small  variation  of  G(z)  is  best  shown  for  a 

particular  case  by  specifying  values  of  II,  r  ,  ct ,,  and  Lz .  Here  we 

o  d 

wish  to  give  a  more  general  demonstration  based  on  the  Taylor  series 
expansion  of  G  about  Az  ■  0.  From  Eq.  (A-40)  we  obtain: 


1 

.  1 

2  ~ 

r 

i  +  e£nrJ 

o  ' 

mi 

+ - 


(A-41) 


In  the  sequel  we  show  that  r 

o 

Eq.  (A-41)  then  becomes: 


R  . 
o 


For  E  ~  -j  and  y 


l.*; 
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(A-42) 


G(Az) 


.0047 


\z 

The  —  tern  becomes  important  (it  A z  •  20H  *  140  km,  the 


tern  at  Az  a  5r  a  5R 
o  o 


2  4 

jjpj  tern  at  Az  15H  *  105  km,  and  thejpj 

45  Y*  i  km  for  a  sea  level  detonation,  where  is  the  yield  in  megatons, 
•i  n 

Recall  tliat  the  maximum  value  of  jiz|  of  interest  is  about 


a  d  .  /,♦: 

a  V 


o 

Cl, 


f>0  1  ♦  ~  Y*  ^  km,  for  c  ■  .3  kn/scc  and 


T'l  100  sec.  Tims  we  find  that  ror  reasonable  yields  of  more  than 
about  1.6  Mt  at  sea  level  the  second  and  third  teres  in  the  Taylor 
series  expansion  of  G  arc  not  inportant  over  the  range  of  Az  or  interest. 

The  cancellation  of  the  various  factors  which  make  up  G(Az) 
can  also  be  seen  by  examining  Eq .  (A-40)  for  Az  large.  For  example. 


for 


<*T 

Az  large  and  positive  ^ 

n-1  ' 


Az 

211 

a  c  and  the  square  root  of 


-Az 


Az 


this  factor  combines  with  the  factors  e  ***  and  e^11  to  produce  «  very 

H  .036  Az/ll 

weax  exponential  dependence  oi  c  '  ■  c 


The  point  is  not  that  G(y)  in  Eq.  (A-39)  is  exactly  constant. 
Doubtless  there  arc  effects  which  could  be  explored  which  depend  on 
deviations  of  G  from  unity.  The  point  is  that  since  G(y)  is  a  slowly 
varying  function. all  values  of  y  from  zero  to  2,  and  hence  all  values 
of  z  extending  to  many  scale  heights  above  sea  level,  arc  of 
nearly  equal  importance  in  Eq.  (A-39).  Physically  this  means  that  the 
shock  front  extending  from  sea  level  to  many  scale  heights  above  sea 
level  is  of  nearly  equal  importance  in  determining  the  far-ficld  Lamb 
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node  amplitudes.  The  significance  of  this  is  that  for  the  first 
maximum  of  the  far-field  vmveform  at  times  like  t  9  the  peculiar 
features  of  the  shock  near  the  ground  such  as  the  presence  of  a 
precursor  should  not  be  of  particular  Importance. 

Finally,  we  return  to  Kq.  (A- 39)  with  G(y)  "  1: 


Pl. 


_ !!* _  f _ lL_  lp± 

r  'i  / sin  r_’^  T  ^  H  :d2  3  1  Vs 
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-(r  ,5)  t?(r  ,£)  r  ^  A!  +  y  + 

o  +o  o  ^  r*  i  J  , 

»y  .  d 
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_  t 


+  y  (A-4  3) 


The  second  integral  in  Kq.  (A-43)  is  always  unimportant  com¬ 
pared  with  the  first  for  the  following  reasons.  We  have  previously 

Ai 

observed  tliat  for  t  3  i.  the  function  — 

d 

•'y 


•  t 

—  ♦  y  has  a  broad 
d 


plateau  for  0  <  y  s  2.  The  second  integral  is  therefore  approximately 
times  the  first  (assuming  this  is  less  than  unity); 


-  -  r 

_1  O  o 


2c  T 


that  is,  it  Introduces  a  correction  factor  of  about  C 

to  the  far-field  Lamb  mode  amplitude.  According  to  the 

.1/3 


-  /- 2a  2 
/  2  +r  -  r 

i  o  o 


1  + 


2c  T 


results  of  the  next  section  r  >  9|-^ 

o  • 


-  9  e*/3M  Y  1/3  km,  where  Y 

HI  •  I 


is  again  the  yield  in  megatons.  For  H«7  km  we  find  that  — r  ,  which  is 
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aaxialzed  at  z  -  311,  is  always  leas  than  .  Expanding  C 

?  ym 

2  '  -2 
1  +  * 


for  — |  <1  gives  C. 
o 


4roCTd 


which,  for  ct^  *  30  Its,  has  a 


euxiaia  value  of  1  + 


.06 


•  1H 

i/j  •*  *  •  r  r  10  to- 


■m 


Thus  neglecting  the  second  Integral  Eq.  (A-43)  becoaes: 
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^  I  2 n  ^  lP» 


S  sin  r  ■*  '  Y  Ht. 
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O  ▼  O  I 
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o  j  r  1 
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(A-44) 


Note  that  If  we  deternlne  rQ  In  Eq.  (A-44)  In  the  "conventional" 


way,  that  Is  for  a  given  i?  take  t.  (r  ,z)  a/y-l  and  r  a  ~l 

+  0  1M  0  IM 

1/y  -  1 

obtain  PL  a  V  .  This  result  has  been  previously  obtained,  as 

has  the  dependence  P.  a  - * — v  t  ~2^2  by  Pierce  and  Posey ^  and 

'sinr.p  d 
r 


*1/3 


by  Pierce,  Posey,  and  Uiff 


ej 

(4) 


In  the  sequel  we  show  that  the  above  aethod  of  choosing  rQ, 
while  probably  correct  for  sasll  yields,  Is  probsbly  not  correct  for 
large  ones.  The  aaln  effect  of  this  finding  Is  to  introduce  a  slightly 
different  dependence  of  P^  on  Y  for  large  yields. 


A.  7  Determination  of  the  Matching  Radius,  rQ 

We  need  a  physical  criterion  to  obtain  the  catching  radius 
r  which  occurs  in  Eq.  (A-44).  The  basic  concept  in  the  analysis  to 
deteralnc  this  matching  radius  is  as  follows: 

The  low  frequency  components  of  the  near-field  disturbance, 
which  eventually  determine  the  far-ficld  amplitudes,  must  become 
spat tally  separated  froa  the  higher  frequency  components  before  the 
linear  Lamb  mode  propagation  theory  is  appropriate.  While  the  high 
and  low  frequency  components  arc  located  in  the  same  region  of  space 
they  interact  nonlinearly.  In  particular,  the  low  frequencies  gain 
energy  from  the  higher  frequencies  as  is  demonstrated,  for  example,  by 
the  lengthening  of  positive  phase  with  distance. 

One  way  this  spatial  separation  of  frequencies  could  occur 
is  the  following.  As  the  shock  or  Mach  stem  propagates  along  the 
ground  it  will  lift  up  because  of  refraction  associated  with  the 
normal  atmospheric  sound  speed  profile.  This  effect  is  responsible 
for  the  so-called  zone  of  silence  phenomenon  at  intermediate  distances 
from  an  explosion,  if  the  low  frequencies  arc  thought  to  exist  in  n 
Lamb  mode  at  this  tine  then  they  will  continue  to  propagate  along  the 
ground  and  hence  will  become  spatially  separate  from  the  higher 
frequencies. 

Refraction  of  the  weak  shock  nay  be  an  inporta  it  part  of 
source  modeling  but  we  do  not  believe  the  above  description  is  the  key 
to  understanding  how  the  Lanb  mode  is  set  up.  The  difficulties  with 
the  above  explanation  are  first  that  it  begs  the  question,  and,  second, 
that  it  relates  to  an  effect  peculiar  to  the  shock  front  near  the 
ground.  In  view  of  the  fact  tliat  the  shock  extending  up  to  at  least 
about  60  km  altitude  is  effective  in  exciting  the  Lanb  mode  we  doubt 
that  this  type  of  refraction  could  by  itself  achieve  the  spatial 
separation  of  frequencies  required. 

The  nodcl  we  propose,  which  is  discussed  below,  involves  a 
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Hp.it  i  a  1  sop.-ir.it  ion  of  frequencies  within  the  weak  shock  itself.  To 
give  the  argiesent  in  the  siaplest  form  we  assiase  that  all  atmospheric 
properties  arc  Independent  of  range.  Then  according  to  Eqs.  (A— 14 ) 

and  (A-18): 


d 


(A-45) 


According  to  what  has  been  previously  said  this  will  be 
approx lEutely  the  quarter  period  of  the  disturbance  at  range  r. 

If  we  start  with  a  hypothetical  far-ficld  waveform,  and 
extrapolate  hack  toward  the  burst,  the  period  will  decrease  according 
to  hq.  (A-45).  At  the  same  tine  we  nay  calculate  the  shcck  positive 
phase  duration,  which  Is  approximately  the  quarter  period  of  the  near- 
field  disturbance,  at  Increasing  distances  fron  the  burst.  This  is 
done  by  using  Kq.  (A— 36)  applied  in  the  horizontal  direction  for  an 
arbitrary  radius  r: 


t+(r,n)  -  t+(Ro)  '  1  *  e  ln  ?  (A-46) 

^  o.. 

it  turns  out  that  there  are  two  cases: 

(a)  ij  >  t+  for  all  radii  at  which  the  shock  is  weak. 

(b)  x j  <  t+  up  to  rQ,  ln  spite  of  the  fact  that  the  shock  is  weak, 

and  t  >  tA  for  r  >  r  . 
d  -  +  -  o 

Consider  the  latter  case  first.  We  will  show  that  the  radius 
r^  has  the  following  physical  meaning:  Up  to  rQ  the  Fourier  component 
of  period  At j  observed  in  the  far  field  is  spatially  located  well 
behind  the  shock  front.  For  radii  greater  than  tq  this  Fourier  com- 
ponent  lias  moved  up  to  the  front  of  the  shock.  While  the  long  period 
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Fourier  components  are  located  behind  the  shock  front  they  gain  energy 

from  the  shorter  period  components  as  Is  manifested,  for  example,  by 

the  positive  phase  lengthening  given  by  Eq.  (A-46) .  Once  these  long 

period  Fourier  components  move  ahead  to  the  shock  fron:  they  propagate 

Independently  of  the  short  period  components. 

On  the  basis  of  this  physical  picture  we  chessc 

r  :  T.(r  )  ■  tA(r  ),  for  case  (b) .  For  case  (a)  we  choose  the  natch- 
o  do  ▼  o 

lng  radius  such  that  a  continuous  transition  is  made  with  case  (b). 

When  this  is  done  it  turns  out  that  case  (a)  corresponds  to  the  method 

1/3 

of  choosing  r  used  by  Pierce  and  Posey,  that  is,  r  a  Y  .  The 
o  o 

criterion  for  case  (b),  however,  is  not  the  conventional  one  and  this 
introduces  a  far-ficld  amplitude  dependence  which  is  not  simply  pro- 
port lonnl  to  yield. 

We  now  give  a  quantitative  statement  of  these  remarks. 

Again  we  take  a  Classtone  pulse  form 

-  _t/t+ 

w(t)  »  *  e  (1  -  t/t+),  t  >  0  (A-47) 

To  see  where  the  low  frequency  components  are  located  within  the  pulse 
we  take  the  Fourier  transform  of  only  that  part  of  the  pulse  between 

l 

the  front  and  a  distance  —  ■  t  behind  the  front: 

c 

t 

S(u)  J  c  14,1  *(t)  dt 

o 


(A-48) 
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for  large  —  the  last  tern  doolnaces  due  to  the  factor 

•  A 


"t/t+ 

e  ,  i  nd  S 


+ 

ITU 


ir“ 


2  •  This  will  be  essentially  the  value  of 


S  when 


U)  > 


iu  ♦ 


(A-49) 


Under  the  assuaption  u  «  ~  this  reduces  to  u  >  e 


-;/t+  j 


or: 


t 


(A-50) 


This  means  that  the  dominant  contribution  to  the  Fourier 
amplitude  of  period  ~  comes  within  a  distance  of  about  ct+  In 

of  the  pulse  front.  For  example,  we  might  have  ~  *  600  sec  and  t+ 

the  order  of  6  seconds,  then  the  dominant  contribution  comes  within 
about  4ct+  of  the  pulse  front.  Under  no  circumstances  is  it  necessary 
to  go  back  more  than  a  few  ct^. 

The  energy  associated  with  a  low-frequency  component  of 
period  4  ?d(r)  (the  period  observed  in  the  far  field)  propagates  at 
a  group  velocity  given  by  Eq.  (A— 10) : 
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c  -  3k  h 


kk 


a  c 


*  c  1 


JL  „ 

3  2t  ,(r)  "kk 
c  l  d  i 


(A-51) 


Similarly  the  energy  associated  with  the  largest  amplitude 


Fourier  component 


w 


1_ 


propagates  at: 


V 
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i 
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fi  'i-i 


,2 


i  *+i 


'kk 


(A-52) 


Justification  for  using  the  Lamb  mode  dispersion  equation  in  deriving 
Eqs.  (A-51)  and  (A-52)  is  given  below. 

At  rQ  therefore  the  low  frequencies  have  moved  ahead  relative 
to  the  main  pulse  by  an  amount: 


\  ~  o  3  1  , 

d(ro)  c  2  2  ‘kk  ’ 

c 


since  Tj  >>  t+.  Using  Eq.  (A-45),  under  the  assumption  that  the  shock 

is  weak  over  most  of  the  distance  to  r  ,  this  becomes: 

o 


d(ro) 


ct .(r  ) 
a _ o 


(A-53) 


According  to  what  has  been  said  before,  we  must  require 

d(r  )  to  be  a  few  times  ct.  since  this  is  how  far  within  the  pulse 
o  +  r 
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the  frequencies  of  interest  I  w 
criterion: 


are  located.  Thus  we  obtain  the 

2Td  i 


t.(r  )  x.(r  )  (A-54) 

f  o  do 

The  sound  speed  occurring  in  Eqs.  (a-51)  -  (A-53)  is  the 
sound  speed  within  the  pulse  which  presumably  is  larger  than  ambient 
due  to  effects  occurring  at  the  front  of  the  weak  shock. 

In  deriving  the  frequency  dependent  group  velocity  in  Eqs. 
(A-51)  and  (A-52)  we  used  the  Lamb  mode  dispersion  equation,  Eq.  (A-10). 
We  now  Justify  this  procedure.  Equation  (A-10)  cannot  be  exactly 
correct  with!  the  shock  because  of  nonlinear  effects.  That  is,  in 
the  Fourier  decomposition  of  the  weak  shock,  rS(w)  e*^x  rS(w) 

is  a  function  of  space.  The  factor  r  is  included  to  account  for 
purely  geometrical  effects.  According  to  Eq.  (A-48)  with  t  ■  00 


r|s(«,r)| 


— 2  +  w 
C+ 


-  2 

a  nrt.  for  u  <  t . 


a  t. 
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A 


1  + 


8*"f- 
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(A-55) 


where  we  have  used  the  fact  that  nrt+  is  a  constant  for  a  weak  shock ^ 
and  have  used  Eq.  (A-46)  to  write  the  explicit  r  dependent  of  t+. 

Because  of  Eq.  (A-55)  it  is  appropriate  to  define  an  imagi¬ 
nary  component  of  k  which  represents  the  spatial  growth  of  frequencies 
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ki  =  h  ls^»r)rl 


V 


1  +  Mn 


=  —  ,  for  r  =  r 
r  o 

o 


(A-56) 


The  real  part  of  k  is  of  course  k  ~  .  If  k  >>  k.  then  the  proup 

r  r  c  r  1 

velocity  will  be  given  by  using  the  Lamb  mode  dispersion  equation. 

r 

This  criterion  amounts  to: 


-  «  r  , 

w  o 


(A-57) 


or  for  the  dominant  Fourier  components  <*)  -  — 

C+ 


ct+  V 


(A-58) 


which  simply  means  that  the  positive  phase  length  must  be  much  less 
than  the  matching  radius.  This  will  be  the  case  since  r^  is  chosen  to 
be  a  radius  where  the  shock  is  weak  and  hence  is  much  larger  than  a 
fireball  dimension  (ct+  being  the  order  of  a  fireball  dimension). 

For  the  very  low  frequencies  observed  in  the  far-field 
(u)  ^  1/Tj(r))  the  criterion  is: 


ct  (r)  <<  r  , 
a  o 


(A-59) 


which  is  not  automatically  satisfied. 

However,  in  deriving  the  matching  criterion,  Eq.  (A- 54),  we 
only  used  the  fact  that  the  very  low  frequencies  have  a  group  velocity 
which  is  much  nearer  the  local  sound  speed  than  that  of  the  high 
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frequencies.  This  will  still  be  the  case  since  both  kf  and  k^  will  be 

small  for  the  low  frequencies  (as  compared  to  k  =  for  the  high 

r  + 

frequencies)  and  hence  the  term  h^  k  In  the  Lamb  mode  dispersion 
equation  (A-10)  may  still  be  omitted  in  determining  dui/dk  for  the  very 
low  frequencies.  Thus  our  use  of  Eq.  (A-10)  in  deriving  the  matching 
criterion  (A-54)  is  justified. 

From  Eqs.  (A-45)  and  (A-46)  this  basic  criterion  (A-54) 

becomes: 


>  3/2 


3h, 


kk 

4  ro 


1  +  6  tnr 

o 


(A— 60) 


According  to  Sachs  scaling: 


1/3 


W  -  *  I  P. 


1/3 


(A-61) 


where  again  using  the  problem  M  results,  c  s  .9  km,  a  ■  1/3  sec,  and 
as  before  3  *  2/3,  when  Y  is  measured  in  kT  and  in  atmospheres. 

3hkk 

Putting  — r-  •  u  we  obtain  from  Eq.  (A-60) : 

c 


l/3°/2 

°'o-“3,5^1  +  s'nr(r)  j 


(A-62) 


•y 

We  wish  to  solve  Eq.  (A-62)  for  rQ  .  It  is  convenient  to  define 

1/3  *  b 

ro  |  Pfc  \ 

the  parameter  x  -  ~  I  ~  I  80  that  E9*  (A-62)  becomes: 
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(A-63) 


or: 


,  .3/2  3/2 

- * - 

1  u3  /  (1  +  gilnx) 


so  chat: 


r 
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1/2 


3/2 


(1  +  £nnx) 
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(A-64) 


(A-65) 


Y  1 3  l  i 

Now  —  has  a  minimum  value  at  x  ■  exp  -r  -  —  |  a  1  where 

*b  t 


Y_ 

P. 


36 


exp 


1  _  — 

36 
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co 


3/2 


f  o 
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,3/2 


Y* 

Pb 


,  for  6  ■  2/3.  At 


this  value  of  x,r  ■  R  ■  e 
o  o 


Y* 


1/3 


y*  *' -  Y 

.  To  continue  the  solution  r  =r- 

PK  °  PK 

b  |  b 


to  values  of  ^  which  are  smaller  than  can  be  obtained  from  Eq.  (A-64) 
Pb 


we  take  r  ■  c 
o 


1/3 


>1/3 


■  .9  |  —  j  kit  for  p-  <  pjj.  This  corresponds 


Y* 


to  the  aforementioned  case  (a)  for  which  t.  >  t.  at  all  radii  for  which 

d  + 


the  shock  is  weak. 


To  see  what  the  net  effect  of  the  altered  criterion  contained 

in  Eqs.  (A-64)  and  A-65)  is  we  return  to  Eq.  (A-43).  The  far-field 

2 

Lamb  mode  amplitude  is  proportional  to  -n(r  )  t,(r  )  r  .  For 

o+oo 

Y  I Y*  l  -2  Y  2 

p-  1  I  jf*  I  t^ls  quantity  is  equal  to  -  10  jT"  km  sec  ,  when  Y  is  in  kT 

b  V  b  /  b 

Y  Y* 

and  P.  in  atmospheres.  For  ^—  <  p^-  we  may  use  the  fact  that  according 

*b  "  *b 
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Thus 


to  weak  shock  then  ry  Tt+  is  independent  of  r 


(2) 


<'„>  'X>  ro  * 10  It  rTF) 

b  ♦  o 


10-*  (i-  l  *  8  t„  xf 


Y  Y* 

The  correction  factor  to  the  far-field  amplitude  when  r—  >  rr 

b  rb 

therefore  be  defined  for  6  ■  2/3  as: 


*  2  )  ** 

Q  1  +  —  tn  x  J  ,  x  >  1 


where  x  is  given  by  Eq.  (A-64) 


‘Y  l  _  j  Y^ 
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p-  1  pbi 


.1/2 


i  +  4  4n  x) 


Tfu 


Using  Eqs.  (A-66)  -  (A-68)  we  uy  rewrite  Eq.  (A-44) 


(A-66) 

nay 

(A-67) 

(A-68) 

as: 


-  98  - 


2 

where  we  have  used  II  •  c  /yg.  The  corresponding  matching  radius  is 
given  by  Eq.  (A-60): 


r 

o 


c 


1/3 


x3/2 

2C  TyZ 

(1  +  -  inx) 


(A-70) 


and  the  critical  yield  to  pressure  ratio  is: 


3hkkC 

4  3 
c  a 


400 


<hkk> 


3/2 


Mt 


sec 


3/2 


atmosphere  ^9/2 


(A— 7 1 ) 


for  c  ■  .9  km,  c  ■  .34  kta/sec,  and  a  ■  1/3  sec.  Rnsed  on  a  particular 
model  atmosphere  and  direction  of  propagation  chosen  to  represent  the 
signal  received  at  Berkcly,  California,  following  the  llousatonic 
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/  j  ^  IqB 

ilfton.it ion  .it  Johnston  Island  Pierce  and  Posey  find  5  *13 


(Y*  t 

—  I  =  19  — — - - -  .  For  a  different 

P*  |  a  tao sphere 

set  of  atmospheric  conditions  near  the  burst  point  or  a  diff»rent 
direction  of  propagation  h^  and  thus  |  jjj-  j  would  be  different. 


We  define  the  scaled  yield  to  pressure  ratio: 


-  1  L* 

%  !p*b 


Q2-  1 


(A-72) 


where  we  linve  used  i.qs.  (A-67)  and  (A-68)  to  eliainate  x. 
the  scaled  matching  radius  is: 


Siailarly 


(A-73) 


For  Q  <  1.6  which  corresponds  to  T  <  4  we  may  derive  a 
simple  analytic  form  for  Q(T)  from  Eq.  (A-72): 


in  f  •  |  A  “■  -1  -  in  Q 


|  ^f1  '  N*l> 


2  j J 


i  (Q-D2 


(A-74) 


and  therefore: 


q  ■>  l  ♦ 


2»n  r 


(A-75) 
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which  gives  the  correction  factor  directly  in  terras  of  the  scaled 
yield  to  pressure  ratio. 


(A-73). 


r 

In  Fig.  A-4  we  plot  Q  and  —  vs.  T  based  on  Eqs.  (A-72)  and 

o 


In  ord  to  utilize  Eq.  (A-75)  or  Fig.  A-4  for  a  particular 


detonation  it  is  necessary  to  calculate  the  critical  yield  to  pressure 
/  Y*  1 

ratio  |  jp*  .  Tills  quantity  involves  the  dlspcrsiu  coefficient  h^ 
near  the  detonation  which  can  only  be  obtained  from  numerical  calcu¬ 


lations  based  on  a  detailed  atmospheric  model.  It  is  important  to 
note  that  I  tit  I  will  differ  in  different  directions  from  the  burst. 
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APPENDIX  B:  VARIATION  OF  RAYLEIGH  WAVE  AMPLITUDE  WITH  YIELD 
AND  HEIGHT  OF  BURST  FOR  INTERMEDIATE  ALTITUDE 


NUCLEAR  DETONATIONS 

B.l  Introduction 

The  shock  wave  from  an  atmospheric  explosion  can  excite 
Rayleigh  waves  when  it  strikes  the  earth.  In  this  appendix  we  derive 
a  theory  viiich  predicts  shock  properties  (overpressure  and  positive 
phase  length)  on  the  ground  beneath  a  nuclear  detonation.  Using  this 
theory  we  then  compute  the  source  strength  for  Rayleigh  wave  excitation 
as  a  function  of  yield  and  height  of  burst  for  intermediate  altitude 
detonations.  A  more  precise  description  of  the  altitude  regime  which 
we  treat  is  given  in  Table  B-2  and  "source  strength"  is  defined  by 
Eq.  ( B  —1 7 )  of  the  sequel.  We  do  not  consider  any  details  of  Rayleigh 
wave  propagation,  so  that  we  are  unable  to  provide  absolute  values  of 
far-field  Rayleigh  wave  amplitudes  or  phases.  The  present  analysis  is 
directly  applicable  only  to  estimating  relative  yields  or  heights  of 
burst  for  situations  in  which  source  and  receiver  locations  are  the  same 
for  two  or  more  detonations.  However,  the  shock  parameters  at  the 
earth's  surface,  which  arc  calculated  in  the  theory,  arc  suitable  as 
Inputs  to  existing  Rayleigh  wave  propagation  models  such  as  those  of 
Harkrlder  and  Flinn^  and  Nickel  and  Whitaker. ^ 

It  is  in  the  inclusion  of  the  influence  of  atmospheric 
structure  on  the  shock  parameters  resulting  at  the  earth's  surface  that 
the  present  treatment  differs  from  previous  work.  Several  authors  have 
utilized  a  homogeneous  model  atmosphere  in  computing  the  time  dependent 
overpressure  at  the  earth's  surface.  The  assumption  of  atmo¬ 

spheric  homogeneity  can  only  be  valid  for  detonation  altitudes 
sufficiently  less  than  an  atmospheric  scale  height. 

More  recently,  Harkrider  and  Flinn^  have  used  a  layered 

model  atmosphere  and  have  assumed  that  a  given  weak  shock  overpressure 

Ap  relative  to  the  ambient  pressure  at  the  source  altitude  p  occurs 

-1/3  c 

at  a  distance  r  which  scales  as  p  .  The  lengtli  of  the  positive 
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phase  of  the  shock  I.  at  the  distance  r  is  also  assumed  to  scale  as  p 

(1)  C 

The  scaling  relations  used  by  Harkrider  and  Flinn  are  usually 
referred  to  as  ordinary  Sachs  scaling. ^  This  scaling  law  for  the 
shock  overpressure  is  known  to  be  invalid  for  weak  shocks  which  propa¬ 
gate  downward  through  several  atmospheric  scale  heights.  Lutzky  and 

Lehto^  have  shown  that  so-called  modified  Sachs  scaling  in  which  r 
-1/3 

scales  as  p  (where  p  is  the  ambient  pressure  at  a  distance  r  from 

the  explosion)  for  a  given  value  of  4p/p  is  a  much  more  accurate  scaling 

law  for  the  overpressure  of  downward  propagating  shock  than  is  ordinary 

Sachs  scaling.  It  should  be  pointed  out  that  Harkrider  and  Flinn^  do 

not  treat  explosions  above  4.88  km  altitude.  For  low  enough  explosion 

altitudes  quantitatively  different  results  due  to  using  different 

ambient  pressure  scaling  laws  should  be  small. 

However,  for  intermediate  altitude  detonations  the  ordinary 

Sachs  scaling  laws  arc  definitely  inadequate  for  predicting  shock 

parameters  at  tne  earth's  surface.  Nor  does  it  appear  that  modified 

Sachs  scaling  is  adequate  for  this  purpose.  This  is  for  two  reasons: 

(1)  The  validity  of  modified  Sachs  scaling  at  relative  overpressures 

much  lower  than  those  which  occur  in  the  calculations  of  Lutzky  and 

Lohto^  is  uncertain.  (2)  Modified  Sachs  scaling  applies  only  to  the 

relative  overpressure  and  not  to  the  positive  phase  length.  In  general 

the  positive  phase  length  is  required  in  order  to  calculate  the  far- 

field  Rayleigh  wave  amplitude,  although  for  detonations  well  above 

(2) 

100  km  altitude  this  may  not  be  the  case. 

In  th'v  absence  of  reliable  height  of  burst  scaling  laws  we 
have  used  the  nunerlcal  calculations  of  Lutzky  and  Lehto^  to  predict 
the  shock  properties  at  the  earth's  surface.  These  authors  have  made 
calculations  of  the  downward -going  shock  from  a  point  explosion  in  an 
ideal  gas,  spherically  symnetric,  exponential  atmosphere.  Several 
factors  must  be  assessed  in  applying  these  idealized  calculations  to 
the  case  of  a  nuclear  explosion  in  the  real  atmosphere: 
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a.  The  assumption  of  an  exponential  atmosphere  limits  the  theory, 
in  its  present  form,  to  the  0-100  km  altitude  interval  where  the 
scale  height  may  be  assumed  constant. 

b.  Application  of  results  obtained  for  a  spherically  symmetric 
atmosphere  (i.e.  one-dimensional)  to  the  real,  horizontally 
stratified  atmosphere  involves  the  assumption  of  independent 
propagation  of  different  portions  of  the  shock  front.  It  has 
been  found  that  differences  between  one  and  two  dimensional 
blast  wave  calculations  in  an  isothermal  atmosphere  are  in 
fact  small.  ' 

c.  The  calculations  of  Lutzky  and  Lehto  are  for  an  ideal  gas  with 
y*1.4.  Near  the  burst  point  this  description  of  the  equation 

of  state  of  air  is  obviously  lacking  because  of  the  high  tempera¬ 
tures  and  pressures  involved.  Numerical  calculations  indicate 
that  the  net  effect  of  this  complication  is  the  following:  a 
nuclear  explosion  at  sea  level  in  real  air  is  0.7  times  as 
effective  in  producing  a  given  weak  shock  overpressure  at  large 
distances  as  the  idealized  point  source  with  y*1»4.  That  is, 
a  0. 7  Mt  idealized  point  source  in  the  y*1>4  medium  produces 

the  same  given  weak  shock  overpressure  as  a  1  Mt  actual  nuclear 
(9) 

burst  in  real  air. 

Even  allowing  these  assumptions,  the  numerical  calculations  of  Lutzky 
and  Lehto  are  not  entirely  sufficient  for  our  purposes  because: 

(1)  The  do  not  extend  to  low  enough  relative  overpressures.  (2)  They 
do  not  explicitly  calculate  the  length  of  the  positive  phase. 

In  the  next  section  we  match  these  numerical  calculations  to 
an  analytic  weak  shock  theory  in  order  to  extend  the  relative  over¬ 
pressure  range  which  may  be  treated.  Also  in  the  next  section  we  use 
this  analytic  weak  shock  theory  to  relate  the  positive  phase  length  to 
the  spatial  derivative  of  peak  overpressure,  a  quantity  which  is  calcu¬ 
lated  numerically.  Thus  by  using  weak  shock  theory  in  conjunction  with 
the  numerical  calculations  we  are  able  to  calculate  shock  parameters  at 
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the  earth's  surface  without  using  Sachs  scaling  laws.  Having  the  shock 
parameters  at  the  earth's  surface  we  then  compute  the  variation  of  far- 
field  Rayleigh  wave  amplitude  with  yield  and  height  of  burst. 

B. 2  Determination  of  Shock  Properties  at  the  Barth's  Surface 

The  calculations  of  Lutzky  and  Lehto^^  are  parameterized  by 
the  quantities  =  r(pc/E)^^  and  ■  h(pc/E)^\  where  E  is  the 
energy  released  by  the  point  explosion,  p^  is  the  ambient  pressure  at 
burst  altitude,  h  is  the  scale  height  and  r  is  the  distance  from  the 
explosion.  Note  that  and  are  scale  height  and  radius  in  units 

1  io 

of  the  length  (e/ p^)  characterizing  the  detonation. 

We  will  match  the  numerical  solution  to  an  analytic  one  at  a 
relative  overpressure  Ap/p  =0.1  (Ap  is  the  overpressure  and  p  is  the 
ambient  pressure  at  a^).  Assume  that  in  the  neighborhood  of  the  match¬ 
ing  point  we  may  write  Ap/p  =  a  o^  ,  where  a  and  n  are  constants.  It 
turns  out  that  matching  of  solutions  can  be  conveniently  done  in  terms 

of  two  quantities: 

This  amounts  to  prescribing  Ap  and  its  first  derivative  at  o^. 

Values  of  o^.n  and  the  useful  quantity  x  »  5^/20^  obtained 
from  the  calculations  of  Lutzky  and  Lehto^^  are  shown  in  Table  B-l. 


TABLE  B-l 

MATCHING  PARAMETERS  OBTAINED 
FROM  NUMERICAL  CALCULATIONS 


°h 

a 

r 

n 

X 

0.05 

0.311 

-4.26 

3.11 

0.10 

0.511 

-3.84 

2.55 

0.20 

0.786 

-3.29 

1.96 

0.50 

1.26 

-2.82 

1.26 

2.0 

2.10 

-1.90 

0.502 

ar  at  Ap/p 


0.1,  and  n 


3  U-  In  ^ 
r  \3or  p  / 
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In  Fig.  B-l,  n  and  are  plotted  vs  x.  We  have  also  included  the 

point  n(x=0)  corresponding  to  a  homogeneous  atmosphere  (a.  “  «)  based 

(9)  1 

on  Figure  3  of  Lehto  and  Larson. 

Now  according  to  the  weak  shock  theory  of  Reed^^ 


L 

L 


r 

o 

L 

o 


I'1 

dr 


(B-l) 


(B-2) 


where  ApQ/po  is  the  relative  overpressure  and  Lq  is  the  length  of  the 
positive  phase  at  <  r  and  where  y  is  the  ratio  of  specific  heats. 

From  Eq.  (B-l)  assuming  an  exponential  atmosphere 


(B-3) 


where  h  is  the  atmospheric  scale  height. 
Using  Eq.  (B-2)  we  find 


Xti  1°  L_ 

-s  i  p  r  \l  p  L 

3  ,  L  ro  V  r  o 

77  ln  r  - - ~2 - 

°  L 

L 


y+1  Ap  .1 
4y  p  L 


where  we  have  used  (B-l)  to  eliminated  Ap  /p  . 

o  o 


( B— 4 ) 
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FIGURE  B-l  SCALED  SCALE  HEIGHT  crh  AND  HATCHING  PAPAWETER 
AS  FUNCTIONS  OF  THE  SCALED  MATCHING  RADIUS  X 


Substituting  (B-3)  In  (B-A)  yields: 


_ Xti.  i 

3r  p  r  2h  Ay  p  L 


or 


(B-5) 


,-<•1  Ap 
^  Y  P 


}+lh  + 


rJ 
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In 


A£l 


(B-6) 


Taking  y«l.A  we  find  that  at  Ap/p  ■  0.1  (o^,  ■  o  ): 


L(or) 


r  .  -0-0857  h  x 
1  +  x  +  n 


(  B— 7) 


While  from  (B-6) 


Y+l  ap  _r 

P  Lo 


-2  1  +  x  + 


U-8) 


so  that  from  (B-2) 


L 
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1-2  (1  +  x  +  n) 


-v  f  £(r )  dr 
\  p  r 


|^1  +  x  +  n  oX 


r 


LE!(X)  " 


(B— 9) 
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where  i  is  the  exponential  integral.  Using  (B-7),  the  length  of  the 
positive  phase  at  distance  r  is  found  to  be 


L 


-0. 083  7  li  x 
1  +  x  +  n 


'  *  x  V  ' 

l-2i  +  x  +  n  e  E,(x)  -  C.  ~r 
to  ^  1  X  /h 


r>r 


(B-10) 


Ktnnllv,  substituting  (B-9)  into  (B  — 1 )  and  taking  Ap  /p  "0.1 

o  o 

we  find 


An. 

P 


p 

0.2h  X-  ex  ~  — 

r  .  p 


\  1-2  1  +  x  +  n  e 
^  * 


Vx)  -  h{k 


(B-ll) 


K  >r  die  class  of  detonations  which  we  treat,  the  quantity 
L^(r/2h)  can  be  neglected  in  (B-10)  and  (B-ll).  This  is  a  valid 
approx imat ion  when  the  relative  overpressure  at  the  ground  is  somewhat 
smaller  tlian  0.1  (corresponding  to  a  reflected  overpressure  at  the 
ground  of  about  200  mb).  In  terms  of  Eq.  (B-10)  this  approximation  is 
related  to  the  fact  that  the  length  of  the  positive  phase  becomes 
"fro7.cn"  for  downward  propagation  at  a  value  independent  of  distance 
from  the  burst  point.  It  is  convenient  for  later  use  to  define  the 
dimensionless  functions 


f(x)  = 


0.2  x  e 


'  1  -2  [l  +  x  +  n  j  eX  E^x)  - 


Apr 


PCP  h 


«(x) 


-0.0857  x  f  f 
1  +  x  +  n  \  L 


1  +  x  +  n  ^  e"  E 
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L 
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(B— 1 2) 


(B-13) 
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The  functions  f(x)  and  g(x)  are  shown  in  Figs.  H-2  and  B-3. 
Lquations  (B-12)  and  (B-13)  fona  the  basis  for  the  subsequent  analysis 
in  which  we  consider  the  variation  of  Rayleigh  wave  amplitude  with 
yield  and  height  of  burst. 

First,  however,  in  concluding  this  section  it  is  useful  to 
summarize  the  various  factors  which  determine  the  altitude  regime  to 
which  the  present  theory  is  applicable.  According  to  what  has  been 
said  previously,  the  point  source  equivalent  energy  E  of  a  sea  level 
nuclear  explosion  of  yield  Y  is  0.7  Y.  We  therefore  redefine  the 
parameter 


(B-14) 


Our  treatment  is  essentially  limited  to  the  range  of  values 
of  used  in  the  numerical  calculations  of  Lutzky  and  Lehto^ 
(0.05-2.0).  For  a  given  yield,  Eq.  (B-14)  therefore  places  constraints 
on  the  relevant  pressure  or  altitude  interval: 


0.685  Y(kT)  >  p(mb)  >  1.07  x  10**  Y(kT) 


(8-15) 


There  are  two  additional  restrictions  on  the  relevant  alti¬ 
tude  interval  which  have  been  previously  mentioned:  (1)  The  theory  is 
limited  to  detonations  below  100  km,  because  above  this  altitude  the 
isothermal  atmosphere  approximation  is  invalid.  (2)  We  are  also  limited 
to  sufficiently  high  altitude  detonations  that  the  relative  overpressure 
on  the  ground  is  smaller  than  0.1 

Putting  all  of  these  factors  together  for  the  values  at 
our  disposal,  we  construct  Table  B-2,  which  shows  the  altitude  interval 
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FIGURE  B-2 :  THE  FUNCTION 


figure  h-i:  ’he  Fur.cno*. 
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0  ^ 


wltl cli  can  be  treated  for  several  yields.  In  deriving  Table  B-2,  we  have 
used  the  1962  Standard  Atmosphere  and  assumed  a  scale  height  h*7  km. 


tabu::  b-2 

ALTITUDE  INTERVALS  TO 


WHICH 

THE  THEORY  APPLIES 

Yield 

Altitude  Interval 

10  kT 

33.8  -  93.5  km 

100  kT 

18.6  -  92.5  km 

1  Mt 

30.7  -  80.0  km 

H. 3  Source  Strength  for  Rayleigh  Wave  Excitation  as  a  Function  of 
Yield  and  Height  of  Burst 

(4) 

Toksoz  and  Ben-Menahem'  have  considered  the  far-fleld 
Rayleigh  wave  amplitudes  excited  by  a  point  source  in  a  homogeneous 
atmosphere.  For  three-dime  slonal  propagation  in  the  absence  of  r  n- 
1 inear  effects,  the  source  parameters  contribute  to  the  amplitude 
almost  entirely  through  a  factor 


1  L(u>)  | 


CD 

J  6p(t)  e~lut  dt 

aWD 


(B-16) 


which  is  the  Fourier  transform  of  the  point  source  pressure  fluctuation 
5p(t).  As  described  by  the  above  authors  the  far-field  Rayleigh  wave 
amplitude  is  proportional  to  |L(u)|  and  also  to  a  factor  which  lies 
between  unity  and  exp  -k^  r,  where  k^  is  the  imaginary  component  of  the 
seismic  disturbance  wave  number,  and  r  is  the  height  of  the  source  above 
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ground.  These  are  the  only  factors  which  depend  on  source  paramet'-rs. 

For  the  atmosphere-lithosphere  contrast  k^  ~  (2*/a^)  x  10  \  where 

is  the  seismic  wavelength.  Thus,  exp  -k^  r  is  essentially  unity,  even 

when  r  is  comparable  to  A^,  and  we  may  simply  take  the  Rayleigh  wave 

amplitude  to  be  proportional  to  L(o)  with  no  other  factors  which 

depend  on  the  source  amplitude,  time  dependence,  or  height  above  ground. 

In  the  case  of  linear  pressure  pulse  propagation, 

6p(t)  a  6p  (t)r,  where  6p  (t)  is  the  time  dependent  pressure  pulse  on 
8  8 

the  ground  a  distance  r  beneath  the  source.  Even  in  the  case  where 
the  pressure  propagation  is  nonlinear,  because  of  pulse  lengthening  and 
dissipative  effects  we  expect  that  the  source  strength  for  Rayleigh 

wave  excitation  will  be  proportional  to  the  Fourier  transform  of  Ap  (t)r. 

8 

The  reason  is  simply  that  the  source  strength  can  only  depend  on 

parameters  at  the  ground  and  not  on  the  previous  history  of  the  pulse 

propagation.  That  is,  we  may  substitute  an  equivalent  linear  point 

source  pressure  fluctuation  proportional  to  6p  (t)r  for  the  actual 

8 

nonlinear  disturbance. 

For  reasons  which  will  soon  be  evident,  we  therefore  define 
the  dimensionless  source  strength  to  be 


S(T) 


2irt 

T  p  h 
8 


on 

J  6p(t)  e  dt 

o 


(B-17) 


where  T  ■  2tt/w  is  the  Rayleigh  wave  period  of  interest,  p  is  the  sea 

8 

level  pressure,  and  h  is  the  atmospheric  scale  height.  Notice  that 
S(T)  does  not  contain  any  dependence  on  source  parameters  except  through 
|L(u) | .  The  far-field  Rayleigh  wave  amplitude  is  therefore  proportional 
to  S(T)  with  no  other  factors  which  depend  on  source  amplitude,  time 
dependence,  or  height  above  ground. 
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To  continue  we  assume  a  pressure  pulse  at  the  ground 


V‘> 


Ap  (l 
\ 


(B-18) 


where  Ap  is  the  peak  overpressure  and  t  is  the  positive  phase  duration. 
The  weak  shock  theory  employed  in  the  previous  section  should  be  Inde¬ 
pendent  of  pulse  shape  since  nonlinear  effects  arc  governed  by  what 
happens  at  the  shock  front  and  not  the  shock  interior.  Thus,  Eq.  (B-17) 

becomes 


S(T) 


Ap  r 


.2 


2*t 


P.h 


+  - 


(B-19) 


Now  t+  *  L/c,  where  c  is  the  sea  level  sound  speed  and  L,  as 
lr  the  previous  section,  is  the  positive  phase  length.  Thus,  in  terms 
of  the  functions  f(x)  and  g(x)  of  the  previous  section,  Eq.  (B-19)  can 
be  written 


for  h*7  km,  and  c*0.34  km/sec.  In  writing  (B-20)  we  have  explicitly 
noted  that  the  ri^ht  sides  of  (3-12)  and  (B-13),  defining  f(x)  and 
g(x),  are  to  be  evaluated  at  the  ground  where  p  ■  p  . 


(ID 
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According  to  Eq.  (B— 14) 


3/2 

'iiii  /  ! Ji)  y1/2  -  9.1975  x  10-3  c3/2  Y1/2 

V  Pg  )  •  h 


(B-21) 


for  p  *  1013.25  mb.  Thus,  (B— 20)  becomes: 

9.1975  x  l(f3  o3/2  Y1/2  f (x) 

S  -  - - - — -  (B-22) 

!  +i _ i _ I 

^ 129.3b  g(x)J 

We  generate  curves  of  source  strength  S  versus  height  of  burst 
r  in  the  following  way:  Picking  Y  and  o.  determines  pc  and  her.ee  r 
through  Eq.  (B-21)  (again  using  the  1962  Standard  Atmosphere).  The 
quantity  x  is  determined  from  Fig.  B-l  or  Table  1J-1,  and  f(x)  and  g(x) 
for  use  in  Eq.  (B-22)  are  found  from  Figs.  B-2  and  B-3. 

Results  for  yields  Y«10  kT,  100  kT,  and  1  MT  and  Rayleigh 
wave  periods  T  ■  10,  20,  and  40  seconds  are  shown  in  Figs.  B-4,  B-5, 
and  B-6.  The  yields  we  refer  to  are  effective  yields  which  at  high 
altitudes  may  be  less  than  the  nominal  yields  due  to  radiative  or  other 
energy  loss  mechanisms.  There  is  a  maximum  in  source  strength  which 
occurs  at  about  80-85  km  for  10  kT,  65-70  km  for  100  kT,  and  50-60  for 
1  MT.  The  strength  of  the  maximum  is  greatest  for  the  larger  Rayleigh 
wave  periods  and  for  the  shorter  periods  the  maximum  may  disappear 
entirely.  The  altitude  at  which  the  maximum  occurs  increases  slightly 
with  Rayleigh  wave  period. 

The  variation  of  source  strength  with  altitude  can  be  appre¬ 
ciable.  For  example,  the  source  strength  for  a  Rayleigh  wave  period  of 
40  seconds  and  a  yield  of  100  kT  increases  by  almost  a  factor  of  3  as 
the  height  of  burst  is  raised  from  20  km  to  70  km. 
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FIGURE  B-A  =  SOURCE  STRENGTH  FOR  RAYLEIGH  WAVE 

EXCITATION  VS  HEIGHT  OF  BURST,  Y=IOKT 


FIGURE  B-5  SOURCE  STRENGTH  FOR  RAYLEIGH  WAVE  EXCITATION 
T  =  10  sec.  VS  HEIGHT  OF  BURST,  Y  =  100  KT 
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Figures  B-7  and  B-8  show  source  strength  plotted  against  yield 

for  Rayleigh  wave  periods  of  10,  20,  and  40  seconds  and  burst  heights 

of  40  and  80  km.  Within  the  range  of  parameters  treated,  the  source 

1/2  2/3 

strength  increases  with  yield  at  a  rate  between  Y  and  Y 

B.4  Conclusions 

a.  Rayleigh  wave  amplitudes  from  intermediate  altitude  detonations 
vary  with  height  of  burst  as  well  as  yield. 

b.  The  variation  with  yield  appears  to  lie  between  the  1/2  and  the 
2/3  power,  depending  on  burst  height  and  the  Rayleigh  wave  period 
being  considered. 

c.  For  a  given  yield  the  Rayleigh  wave  amplitude,  at  least  for 
periods  of  20  seconds  or  more,  has  a  maximum  at  a  particular 
burst  height.  The  maximum  occurs  at  80-85  km  for  10  kT,  65-70  km 
for  100  kT,  and  50-60  km  for  1  Mt.  The  strength  of  the  maximum 
is  greatest  for  the  longer  Rayleigh  wave  periods  and  the  altitude 
at  which  the  maximum  occurs  increases  slightly  with  increasing 
Rayleigh  wave  period. 
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GLOSSARY  OF  SYMBOLS  USED 


Section  2  and  Appendix  A 

A  amplitude  function  of  Lamb  mode 

A^  Airy  function 

A^  derivative  of  the  Airy  function 

c  sound  speed 

C  correction  factor  due  to  Az  <  0  integration  in  Eq.  (A-43) 

d  distance  moved  by  low  frequencies  relative  to  high  frequencies 

F  function  containing  time  and  horizontal  distance  dependence 

of  Lamb  mode 

G  function  defined  by  Eq.  (A-40) 

h^  dispersion  coefficient  due  to  deviations  of  wind  and  sound 
speeds  from  vertically  averaged  values 

H  scale  height 

k  wave  number 

imaginary  part  of  k 

kf  real  part  of  k 

£.  distance  behind  shock  front 

ambient  pressure  at  burst  altitude 

P£  critical  ambient  pressure  for  breakdown  of  yield-amplitude 

proportionality 

P^  Lamb  mode  pressure  amplitude 

Pg  sea  level  pressure 

pressure  amplitude,  first  peak  to  trough 
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AP  overpressure 

Q  large  yield  amplitude  correction  factor 

r  great  circle  distance 

r  earth  radius 

e 

rQ  cylindrical  radius  for  matching  weak  shock  and  Lamb  mode 

r  parameter  defined  by  Eq.  (A-14) 

R  distance  from  burst 

Re  "real  part  of" 

R  small  radius  where  ir(R  )  -  .1  and  atmospheric  inhomogeneity 

is  unimportant 

S (to)  Fourier  transform  of  Glasstone  pulse 

t  time  after  signal  arrival 

t+  positive  phase  duration 

t  parameter  defined  by  Eq.  (A-13) 

T  period  of  far-field  disturbance 

V  group  velocity 

6 


Y  yield 

yield  in  megatons 

Y*  critical  yield  for  breakdown  of  yield-amplitude  proportionality 

z  altitude 
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z 

m 

z 

Az 


a 


0 


Y 

r 


altitude  of  rigid  ceiling  in  modal  analysis 
height  of  burst 

vertical  distance  relative  to  height  of  burst 
lowest  altitude  to  which  shock  extends 
highest  altitude  to  which  shock  extends 

t+  at  it  =  .1  for  1  kT  at  sea  level 

..  rr(R  )  R 
Y+l  o  o  =  o 

2y  ct+(Ro)  "  P 
ratio  of  specific  heats 

scaled  yield  to  pressure  ratio  defined  by  Eq.  (A-72) 


c 


e  radius  at  which  tt  *  .1  fcr  1  kT  at  sea  level 

Tr(r  ,z,t  )  relative  overpressure  at  r  ,z,t 
o  o  o  o 

tt  peak  relative  overpressure 

t  time  after  passage  of  shock  front 

t  waveform  arrival  time 

a 

characteristic  dispersion  time 


waveform  function  of  Lamb  mode 

A 

^  Fourier  transform  of  ij/(r) 

4*  sum  of  modes  other  than  the  Lamb  mode 

a)  frequency  of  disturbance 
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Section  3 


—  Bernoulli  numbers 
c  sound  speed 

d  distance  between  ground  reflections 

exponential  integral  given  by  Eq.  (3-14) 

E^  exponential  integral  given  by  Eq.  (3-13) 

f  function  giving  relative  overpressure  vs.  distance  for 

1  kT  at  sea  level 


H 

H 

I 


i 

m 

n 

nh 

P 

Pb 

AP 


r 

r 

o 

R 


T 


effective  scale  height 
atmospheric  scale  height 
integral  defined  by  Eq.  (3-8) 

contribution  to  I  from  £th  transit  between  ground  and  upper 
boundary 

index  for  each  transit  between  ground  and  upper  boundary 


number  of  reflections  at  upper  boundary 


AP 

P(r) 


value  of  n  for  a  homogeneous  atmosphere 


ambient  pressure 


ambient  pressure  at  burst  point 


overpressure 

radial  distance  from  burst 

value  of  r  at  =  .2 
P 

range 

positive  phase  duration 
far-field  disturbance  period 
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Y  yield 

zq  altitude  of  upper  reflecting  boundary 

z  height  of  burst 

Y  ratio  of  specific  heats 

T  Euler's  constant 


Section  4  and  Appendix  B 


a 

c 

E 

E1 

f(x) 

g(x) 

h 

L(u>) 

L 


constant  relating  relative  overpressure  to  radius  at 


Ael 


.1 


Pb’Pc 


g 


‘  o 

Ap 

AP, 


sound  speed 

point  source  energy 

exponential  integral 

dimensionless  function  defined  by  Eq.  (4-8) 
dimensionless  function  defined  by  Eq.  (4-7) 
scale  height 

Fourier  transform  of  point  source  pressure  fluctuation 
positive  phase  length  at  r 

o 

positive  phase  length  at  ^  ■  .1 
starting  value  defined  by  Eq.  (4-10) 
ambient  pressure  at  r 
ambient  pressure  at  burst 
sea  level  ambient  pressure 
ambient  pressure  at  rQ 
peak  overpressure  at  the  ground 
overpressure  at  r 
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point  source  pressure  fluctuation 

pressure  fluctuation  at  the  ground  from  a  print  source 

radius  from  burst  point,  also  height  of  burst 

arbitrary  radius  loss  than  r,  r  in  particular 

radius  at  ^  ■  .1 
P 

dimensionless  source  strength  defined  by  Eq.  (4-4) 
time 

positive  phase  duration  at  the  ground 
r 

o 

2h 

yield 

ratio  of  specific  heats 

scaled  scale  height  defined  by  Eq.  (4-11) 

scaled  radius 

scaled  radius  at  which  ^  •  .1 

P 

Rayleigh  wave  period 
frequency  of  disturbance 


